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Abstract 

The  physical  and  chemical  mechanisms  responsible  for  the  removal  of  ammonia  from  the  atmosphere  have  been  reviewed.  Capture 
by  atmospheric  moisture  (clouds,  rain,  fog),  surface  water  (rivers,  lakes,  seas),  and  deposition  on  vegetation  and  soil  constitute  the  main 
pathways  for  ammonia  removal  from  the  troposphere.  Ammonia  catalyzes  the  atmospheric  oxidation  of  sulfur  dioxide  to  sulfur  triox¬ 
ide  and  reacts  rapidly  with  acidic  components  of  the  atmosphere  (sulfuric,  nitric,  and  hydrochloric  acids).  The  ammonium  salts  formed 
are  the  main  components  of  smog  aerosols  and  thus  affect  the  opacity  of  the  atmosphere  and  the  earth  radiation  budget.  Slow  oxida¬ 
tion  of  ammonia  in  the  atmosphere  plays  only  a  minor  role  in  its  removal.  The  data  obtained  for  ammonia  reactions  under  normal  atmo¬ 
spheric  conditions  are  generally  applicable  to  model  chemical  reactions  occurring  during  massive  release  of  ammonia  in  the  atmosphere, 
provided  the  impact  of  high  ammonia  concentration  on  the  mass  transfer  processes  that  control  some  of  these  reactions,  are  taken  into 
account. 

©  2004  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Ammonia;  Atmospheric  chemistry;  Plume  dispersion;  Modeling;  Catastrophic  release 


1.  Introduction 

The  long-standing  concern  for  public  safety  in  cases  of 
accidental  releases  of  toxic  materials  from  chemical  plants 
has  gained  additional  momentum  from  recent  fear  of  terror¬ 
ist  attacks  on  chemical  plants  or  the  chemical  transportation 
system.  In  cases  of  accidental  or  malicious  massive  release 
of  hazardous  chemicals  in  the  environment,  it  is  essential 
that  first  responders  have  an  accurate  forecast  of  the  spatial 
and  temporal  distribution  of  released  toxic  chemicals  to  de¬ 
cide  on  population  evacuation  needs  and  to  plan  remedial 
measures. 

A  number  of  mathematical  models  exist  that  adequately 
describe  the  dispersion  of  a  plume  of  chemicals  based 
on  prevailing  meteorological  conditions,  terrain  profiles, 
and  physical  properties  of  released  chemicals.  Recent  ad- 
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vances  in  computer  technology  allow  these  models  to  be 
run  quickly  and  effectively  from  laptop  computers,  thus 
making  them  readily  available  to  civil  authorities  directing 
the  response  to  chemical  incidents.  However,  these  models 
generally  lack  an  accurate  description  of  chemical  trans¬ 
formations  that  the  released  chemicals  will  undergo  in  the 
atmosphere.  As  a  starting  point  in  our  effort  to  correct  this 
deficiency,  we  have  prepared  this  review  of  the  present 
knowledge  of  atmospheric  reactions  of  ammonia,  placing 
particular  emphasis  on  kinetics  of  the  reactions  involved 
and  in  identifying  missing  elements  for  a  comprehensive 
incorporation  of  chemical  atmospheric  behavior  of  am¬ 
monia  as  part  of  a  realistic  plume  dispersion  model.  The 
selection  of  ammonia  stems  from  its  economic  importance 
and  widespread  distribution  throughout  the  country  for  in¬ 
dustrial  and  agricultural  use.  Following  this  critical  review 
of  the  significant  body  of  knowledge  on  the  atmospheric 
chemistry  of  ammonia,  recommendations  will  be  presented 
for  additional  investigations  required  to  generate  the  data 
needed  for  an  adequate  description  of  the  chemical  be- 
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havior  of  high  concentration  plumes  of  ammonia  in  the 
atmosphere. 

2.  Aqueous  chemistry  of  ammonia 

2.7.  Dissociation  of  ammonia  in  aqueous  solutions 

In  aqueous  solutions,  NH3  is  partially  hydrolyzed  to  am¬ 
monium  hydroxide  according  to  Eq.  (1): 

(NH3)aq  +  H20<B.NH4+  +  0H-  (1) 

with  an  equilibrium  constant  K\,  =  1.774  x  10“^  at  25  “C 
(298 K)  [1].  Eq.  (2)  represents  the  best  fit  of  the  values 
of  K\i  determined  by  Bates  and  Pinching  over  the  0-40  °C 
temperature  range  [2]. 

4411  025 

In  =  16.9732 - - 0.04407  (2) 

where  the  temperature,  7,  is  in  Kelvin  (7  (K)  =  273.15  -f 
°C).The  pH  of  the  aqueous  solution  strongly  influences  the 
reactions  establishing  equilibrium  1.  In  solutions  of  ammo¬ 
nia  in  water  or  in  aqueous  solutions  at  neutral  and  alkaline 
pH,  the  following  rate  constants  have  been  measured  using 
electric  field  pulse/conductimetric  techniques  [3,4]: 

(NH3)aq  +  H2O  NH4+  +  OH", 

ifc295K  =  5  X  lO^s-i  (3) 

NH4+  +  OH"  (NH3)aq  +  H2O, 

A:295K  =  3  X  (4) 

In  the  acid  pH  range,  similarly  fast  reaction  rates  were 
determined  using  nuclear  magnetic  resonance  tech¬ 
niques  [5]. 

(NH3)aq  +  H+  NH4+,  kigsK  =  4.3  x  M"*  s"^ 

(5) 

NH4+  ^  (NH3)aq  +  H+,  *295 K  =  24.6  8“*  (6) 


2.2.  Solubility  of  ammonia  in  aqueous  solutions 

2.2,  L  Solubility  in  pure  water 

Equilibrium  partial  pressure  data  for  NH3  above  pure 
aqueous  solutions  are  available  over  a  wide  range  of  tem¬ 
perature  and  concentrations.  A  complete  review  and  critical 
evaluation  of  the  data  was  conducted  in  1989  and  published 
by  Clegg  and  Brimblecombe  [2]. 

(NH3)g  ^  (NH3)aq  (7) 

(NH3)aq  +  H20-e-NH4++0H-  (1) 

The  solubility  of  a  weak  electrolyte  such  as  NH3  can  be 
described  in  two  stages;  as  dissolution  followed  by  base  dis¬ 
sociation,  yielding  the  following  equations  for  thermody¬ 
namic  Henry’s  law  constant  Ku  (molkg“^  atm“^),  and  base 
dissociation  constant  Kb  (mol  kg"  ^): 

(fromEq.(7))  (8) 

PNHs 

^  mum^Youmn 

)^WnH3<2w 

where  m  denotes  molality  (mol  kg"^),  p  the  partial  pressure 
(atm),  Y  denotes  the  activity  coefficient,  and  the  activity 
of  water. 

To  determine  Henry’s  law  constant  experimentally,  two 
quantities  are  measurable:  the  partial  pressure  of  NH3,  pnh3» 
over  the  solution,  and  the  total  concentration  of  ammonium 
species  in  solution:  NH3  H-NHJ  =  //tot-  To  compute  Kh,  the 
activity  of  (NH3)aq  must  be  estimated  from  the  measured 
A^tot  and  Kq,  applying  appropriate  approximations  for  the 
activity  coefficients  of  the  species  involved,  depending  on 
the  characteristics  of  the  aqueous  solution  used.  For  exam¬ 
ple,  if  the  solution  is  dilute  enough,  the  activity  coefficients 
can  be  assumed  to  be  equal  to  1.  Many  of  the  differences 
for  the  value  of  Kh  found  in  the  literature  can  be  explained 
by  the  difference  in  the  choices  made  in  the  estimation  of 
the  activity  coefficients.  The  most  recent  determinations  of 
Kh  are  summarized  in  Table  1 . 

The  values  of  Kh  reported  by  Hales  and  Drewes  [6]  are 
clearly  higher  than  those  of  other  investigators,  a  situation 


Table  1 

Determination  of  Henry’s  law  constant 


Reference 

\nKii  =f(X)  (Kh:  kgatmmor*) 

Ku  (kg  atm  mol  ^), 
273.15  K  (O'^C) 

^Th  (kg  atm  mol  *), 
283.15 K  (10°C) 

Ku  (kg  atm  mol  ^), 
298.15  K  (25  °C) 

Hales  and  Drewes  [6] 

231.92 

144.10 

75,45 

Edwards  et  al.  [7] 

160.559  -  (8.621.06/7) 

217,72 

126.19 

60.80 

-(25.6767  In  7)  +0.0353887 

Dasgupta  and  Dong  [8] 

-9.70  +  (4092/7) 

196,52 

115.78 

55.96 

Clegg  and  Brimblecombe  [2] 

-8.09694  +  (3917.507/7)  - 

218.62 

127.66 

60.72 

(0.003147) 

Shi  et  al.  [9] 

-3.221  +  (1396/7) 

77.58 

51.19 

28.92 

Thermodynamic  value  calculated 

-  9.94  +  (4166/7) 

202.58 

118,28 

56.42 

from  NBS  data  [10] 

Thermod)mamic  value  calculated 

57,64 

from  Wagman  data  [1 1] 

. 
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resulting  most  likely  from  poor  management  of  tempera¬ 
ture  changes  and  control  in  the  experimental  procedures, 
as  discussed  by  Dasgupta  and  Dong  [8],  At  the  other  ex¬ 
treme,  Shi  et  al.  [9]  report  values  that  are  much  lower  than 
other  researchers.  Whereas,  other  investigators  ensure  the 
systems  are  at  equilibrium  before  measuring  ammonia  pres¬ 
sure  above  its  aqueous  solution,  Shi  et  al.  used  a  droplet 
train  flow  reactor  in  which,  from  their  own  admission  ([9], 
p.  8816)  “under  the  conditions  of  the  droplet  experiments, 
equilibrium  is  not  attained  and  the  pH  near  the  surface  is 
time  dependent”.  Their  values  for  are  derived  from  the 
model  of  kinetics  of  the  process  and  from  the  kinetics  data 
that  show  significant  uncertainty,  particularly  at  the  high 
end  of  their  temperature  range.  Moreover,  they  used  a  train 
of  aqueous  sodium  hydroxide  solution  (pH  13)  droplets  to 
maintain  alkaline  pH  during  the  experiment.  The  salting-out 
effect  of  sodium  hydroxide,  reportedly  10%,  was  taken  into 
account  using  an  unspecified  methodology  to  estimate  mag¬ 
nitude.  The  accuracy  of  the  values  for  Kh  is  therefore  open 
to  question.  In  a  later  note  [12],  the  authors  recognized  the 
shortcomings  of  their  analysis  and  confirmed  the  reported 
value  for  was  too  low. 

With  the  exception  of  the  two  investigations  discussed 
above,  the  remainder  of  the  investigation  report  values  for 
are  within  less  than  10%  of  each  other.  Dasgupta  and 
Dong  [8]  used  dilute  solution  ammonia  (0.002-0.1  M)  in 
phosphate  buffers  (pH  6-10).  The  activity  coefficient  for  the 
ammonium  ion  was  calculated  using  the  model  of  Laitinen 
and  Harris  [13]  that  accounts  for  the  effect  of  the  overall 
ionic  strength  of  the  solution  and  the  activity  coefficient  for 
(NH3)aq  was  set  as  1.  Edwards  et  al.  [7]  re-analyzed  the 
data  of  Van  Krevelen  et  al.  [14]  using  a  thermodynamic  ap¬ 
proach  where  the  vapor  phase  coefficient  of  ammonia  over 
the  solution  is  estimated  from  the  second  coefficient  of  the 
virial  equation  of  state  for  pure  ammonia.  Given  the  rela¬ 
tively  high  concentration  of  ammonia  used  and  the  very  low 
value  for  (1.774  x  lO"^  at  25  °C),  they  neglected  ammo¬ 
nia  dissociation  and  used  Eq.  (10)  to  estimate  the  activity 
coefficient  of  aqueous  ammonia. 

In  yK  =  2/Snnwn  (10) 

where  subscript  N  is  thereafter  used  to  designate  NH3,  yN 
is  the  activity  coefficient  for  undissociated  aqueous  ammo¬ 
nia,  )0NN  the  two-body  interaction  coefficient  for  aqueous 
ammonia,  and  mN  is  the  molality  of  ammonia  in  water. 

Equating  the  chemical  potentials  of  ammonia  in  the  aque¬ 
ous  solution  and  in  the  gas  phase  and  neglecting  the  vapor 
pressure  of  water  as  much  smaller  than  the  total  system  pres¬ 
sure,  leads  to  Eq.  (11) 

In  =  In  +  2)5nnw2n  (11) 

where  Is  the  fugacity  coefficient  for  gaseous  ammonia: 
cxp{B^P/RT),  5n  second  coefficient  of  virial  equation  of 
state  for  pure  ammonia,  and  vn  is  the  partial  molar  volume 
of  ammonia  (30.0  at  25  °C). 


The  left  side  of  Eq.  (11)  is  computed  from  experimental 
values  and  plotted  against  wn-  The  intercept  gives  the  value 
of  In  and  the  slope  the  value  of  2^nn- 
Clegg  and  Brimblecombe  [2]  re-analyzed  the  data  of  Chen 
et  al.  [15],  neglecting  the  dissociation  of  the  base  in  the 
relatively  concentrated  ammonia  solutions  in  pure  water 
used  by  Chen  et  al.  [15],  considering  the  low  value  of  ATb 
(1.774  X  10-5  at  25  °C). 

The  values  of  reported  by  these  three  last  investi¬ 
gations  are  in  fairly  close  agreement  with  those  calculated 
using  available  free  energy  and  enthalpy  data  (NBS  data). 
Therefore,  within  a  10%  maximum  accuracy  range,  the  re¬ 
sults  of  these  investigations  are  in  good  agreement  and  con¬ 
sistent  with  those  of  thermodynamic  calculations. 

Considering  that  in  the  intended  application,  high  con¬ 
centration  of  ammonia,  will  be  involved  together  with  water 
with  low  dissolved  salt  content  (rain  drops),  we  recommend 
using  the  relationship  developed  by  Clegg  and  Brimble¬ 
combe  [2]  because  it  was  derived  from  data  obtained  at 
high  concentrations  of  ammonia  in  pure  water. 

2.2.2.  Solubility  in  multi-component  aqueous  solutions 
When  ammonia  is  released  in  the  environment,  it  will 
interact  with  water  sources  containing  significant  quanti¬ 
ties  of  dissolved  chemicals  e.g.,  seawater  or  condensation 
droplets  with  dissolved  atmospheric  pollutants.  The  solu¬ 
bility  of  ammonia  will  be  influenced  by  the  presence  of  the 
dissolved  chemicals.  Two  main  approaches  have  been  used 
to  compute  the  solubility  of  ammonia  in  multi-component 
aqueous  solutions. 

As  discussed  earlier,  Edwards  et  al.  [7]  used  a  thermody¬ 
namic  approach  to  analyze  the  equilibrium  of  ammonia  be¬ 
tween  the  aqueous  solution  and  the  gas  phase  and  derive  a 
general  formulation  for  the  activity  coefficient  of  ammonia 
in  aqueous  solution  represented  by  Eq.  (12): 

where  a  is  the  Debye-Huckel  parameter  (0.391  at  25  °C), 
the  charges  of  the  cation  and  anion  of  compound  I,  I 
the  ionic  strength  of  the  solution  (Emm),  fiij  the  two-body 
interaction  coefficient  of  component  i  with  component y,  and 
mj  is  the  molality  of  component  j. 

This  expression  represents  the  deviation  from  non-ideality 
due  to  specific  interaction  between  the  components  in  the 
solution: 

•  The  first  term  represents  the  contribution  of  electro¬ 
static  interaction  between  charge  ions,  based  on  Debye- 
Huckel’s  theory  of  very  dilute  electrolyte  solutions. 

•  The  second  term  is  a  summation  of  all  the  short-range  two- 
body  interactions  between  the  components  in  the  solution: 
o  molecule-to-molecule  Vander  Waals  interactions; 

o  molecule-ionic  species  interactions  (salting-out  effect); 
o  ions-to-ions  interactions  other  than  electrostatic. 
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Table  2 

Molecule-molecule  interactions  for  single  solute 


Solute,  I 

A’l  (kg  water  mol 

NH3 

0.017 

CO2 

0.010 

SO2 

-0.277 

H2S 

0.005 

Table  3 

Ion-ion  specific  interaction  parameters 

Ion,  I 

(kg  water  mol“^) 

NH4‘*‘ 

-0.12 

HCO3" 

0.25 

COs" 

-0.12 

HS03“ 

0.24 

S03- 

-0.05 

NH2COO” 

0.33 

H+ 

0.10 

OH“ 

0.15 

These  interactions  coefficient,  are  determined  empiri¬ 
cally,  as  illustrated  in  the  preceding  section  for  the  case  of 
ammonia  in  pure  water  solution.  In  this  case,  dissociation  of 
ammonia  was  neglected,  thus  zeroing  the  electrostatic  inter¬ 
action  term,  and  ammonia  was  the  only  species  in  solution, 
thus  fiij  =  and  Eq.  (12)  reduces  to  Eq.  (10).  This  ap¬ 
proach  was  used  to  determine  empirically  the  various  inter¬ 
action  coefficients  reported  in  Tables  2-4. 

The  values  for  the  appropriate  parameters  are  inserted  in 
Eq.  (12)  to  obtain  the  value  for  In  yN  that,  together  with  Eq. 
(13),  is  used  to  characterize  the  distribution  of  NH3  between 
the  gas  phase  and  the  aqueous  solution. 

=KHyHmn  (13) 

where  yn  is  the  molar  fraction  of  NH3  in  the  gas  phase. 

Clegg  and  Brimblecombe  [2]  used  a  like  approach  based 
on  a  more  sophisticated  model  of  Pitzer  [16]  to  calculate 
activity  coefficients  in  a  mixed  electrolyte  solution  using  a 
virial  expansion  of  terms  of  concentration.  The  functions 
involved  are  quite  complex,  particularly  for  ionic  species. 
For  neutral  species  like  undissociated  ammonia,  the  activity 
coefficient  is  given  by  Eq.  (14): 

In  yM  =  2^  +  2^  ^Na"»a  +  2^ 

can 

+  (14) 

i  J 

Table  4 

Molecule-ion  specific  interaction  parameters 


Molecule-ion,  ij  =  nj  (kg  water  mol 


NH3-HCO3- 

-0.031 

NH3-CO3- 

0.068 

NH3-HS03“ 

-0.038 

NH3-S03” 

0.044 

NH3-NH2COO- 

-0.028 

NH3-H+ 

0.015 

NH3-OH- 

0.021 

Table  5 

Interaction  parameters  for  NH3  at  25  °C 


Species,  I 

^N/ 

Mnn/ 

NH3 

0.01472 

Mg?+ 

-0.21 

Ca2+ 

-0.081 

Li+ 

-0.038 

Na+ 

0.0175 

-0.000311 

K+ 

0.0454 

-0.000321 

Nai+ 

0 

-0.00075 

ci- 

0 

OH- 

0.103 

NO3- 

-0.01 

-0.000437 

SOs^- 

0.158 

S04^- 

0.140 

C03^- 

0.180 

0.000625 

N02- 

-0.003 

where  is  the  ammonia-cation  interaction  parameter, 
the  ^Na  ammonia-anion  interaction  parameter,  )0Nn  die 
ammonia-neutral  molecule  interaction  parameter,  /tXNy  the 
ammonia-species  i-species  j  (three-bodies  interaction  pa¬ 
rameter),  and  nti  the  molality  of  species  /. 

Clegg  and  Brimblecombe  [2]  conducted  an  extensive  re¬ 
view  and  analysis  of  the  literature  on  the  solubility  of  ammo¬ 
nia  in  a  variety  of  electrolyte  solutions  and  determined  and 
compiled  values  for  various  interaction  parameters.  Table  5 
summarizes  their  results  for  the  species  of  interest  to  this 
project. 

It  should  be  noted  that  the  values  of  ^nn  found  in  these 
two  investigations  are  in  good  agreement:  0.01472  (Table  5) 
and  0.017  (Table  2).  In  general,  for  solutions  that  are  not  too 
concentrated  (w,-  <  1.5molkg“^),  the  three  bodies  interac¬ 
tion  can  be  neglected. 

Although  more  complex,  Pitzer’s  approach  has  been 
found  repeatedly  to  describe  accurately  the  behavior  of 
complex  electrolyte  solutions  such  as  seawater,  brine,  and 
industrial  solutions.  Its  use  is  recommended  to  estimate  the 
solubility  of  ammonia  in  concentrated  salt  solutions. 

2.J.  Uptake  of  gaseous  ammonia  by  water  surfaces— mass 
accommodation  coefficient 

2,3.  L  Modeling  gas-liquid  interactions 

In  a  recent  paper,  Shi  et  al.  [9]  gave  an  excellent  phe¬ 
nomenological  description  of  the  interactions  between  gas 
and  liquid  with  specific  application  to  the  ammonia-water 
interaction.  The  process  can  essentially  be  described  as  a 
sequence  of  the  following  steps: 

•  The  gas  molecule  is  transported  to  the  liquid  surface,  usu¬ 
ally  by  gas-phase  diffusion. 

•  The  gas  molecule  strikes  the  surface  of  the  liquid  and  is 
thermally  accommodated  as  an  adsorbed  surface  species. 

•  Depending  on  chemical  affinity  of  the  gas  for  the  liquid 
and  on  the  conditions  at  the  surface  of  the  liquid,  the 
following  additional  steps  can  occur: 
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o  the  adsorbed  molecule  can  desorb  from  the  surface  and 
return  to  the  gas  phase; 

o  the  adsorbed  molecule  can  become  dissolved  in  the 
liquid  and  migrate  into  the  bulk  of  the  liquid  by 
liquid-phase  diffusion; 

o  the  adsorbed  molecule  can  undergo  reversible  reaction 
with  the  liquid,  such  as  ionization,  or  irreversible  reac¬ 
tion  with  the  solvent  or  dissolved  species  in  the  solvent. 


adsorbed  molecules  and  to  generate  free  space  at  the  inter¬ 
face  to  accept  new  gas-phase  molecules.  The  correspond¬ 
ing  bulk  uptake  coefficient,  Tb,  takes  into  account  the  ef¬ 
fect  on  gas  uptake  of  Henry’s  law  solubility,  the  chemical 
reaction  rates  in  the  bulk  of  the  liquid,  and  the  diffusion 
coefficients  of  the  species  in  the  liquid.  Mathematical 
description  of  these  phenomena  is  usually  complex  and 
of  limited  practical  utility  for  modeling  purposes. 


The  overall  process  can  be  modeled  using  an  electrical 
circuit  analogue  where  each  step  is  analogous  to  an  elec¬ 
trical  resistance.  The  measured  uptake  of  ammonia,  ymeas> 
is  inversely  proportional  to  the  sum  of  the  resistance  rep¬ 
resenting  each  elementary  step  in  the  process.  To  develop 
a  mathematical  framework  for  the  model,  the  following  pa¬ 
rameters  are  defined: 

•  The  gas  molecule  strikes  the  surface  of  the  liquid  and  is 
thermally  accommodated. 

•  l/Fdiff  represents  the  resistance  to  gas-phase  diffusive 
transfer.  For  droplets,  the  value  is  estimated  using  the 
Fuchs-Sutugin  [17]  formulation: 

1  0.75  +  0.283Kn 

Tdiff”  Kn(l  +  Kn)  ^  ^ 

where  Kn  is  the  Knudsen  factor  defined  as  6D^lcdf,  Dg  the 
gas-phase  diffusion  coefficient  of  ammonia,  df  the  effective 
diameter  of  the  droplet,  and  c  the  mean  thermal  velocity  of 
ammonia  in  the  gas  phase 

The  mass  accommodation  coefficient,  or,  represents  the 
probability  that  a  molecule  striking  the  surface  of  the  liquid 
enters  into  the  bulk  of  the  liquid  phase.  In  the  absence  of  re¬ 
action  at  the  interface,  the  fraction  of  molecules  entering  the 
bulk  of  the  liquid  is  equal  to  the  fraction  of  the  molecules 
striking  the  surface  (the  adsorption  coefficient  5)  minus  the 
fraction  of  molecules  desorbing  from  the  surface.  Consider¬ 
ation  of  the  sorption/desorption  equilibrium  at  the  interface 
leads  to: 

^ _  (1/ ~h  ^des  (16) 

or  Sksdi 

where  k^es  Is  the  rate  constant  for  the  desorption  process,  ^soi 
the  rate  constant  for  the  passage  in  solution  of  the  adsorbed 
molecule,  and  a  is  the  mass  accommodation  coefficient,  is 
thus  a  measure  of  the  ratio  of  the  solution  to  desorption  rate 
constants. 

•  Chemical  reactions  at  the  interface  open  a  new  channel 
for  the  removal  of  the  adsorbed  gas-phase  molecules,  dis¬ 
tinct  from  the  mass  accommodation  coefficient,  a.  These 
reactions  are  accounted  for  by  placing  an  additional  re¬ 
sistance,  l/Fs,  in  parallel.  The  surface  uptake  coefficient, 
Fsy  is  proportional  to  the  ratios  of  the  rate  constant  of  the 
reactions  occurring  at  the  interface  over  the  rate  constant 
of  the  desorption  process. 

•  Solubility  and  reactions  in  the  bulk  of  the  liquid  are 
usually  the  main  avenues  to  deplete  the  interface  from 


The  measured  uptake  of  ammonia,  ymeas*  Is  thus  given  by: 

1  _  J_  1  /  1 

ymeas  ‘fdiff  S  \  -Fs  H"  ""  -H) 

By  applying  semi-empirical  relationships  such  as  the 
Fuchs-Sutugin  equation  (Eq.  (15))  to  estimate  1/Fdifif,  dif¬ 
ferent  pH  conditions,  and  known  reaction  rate  constants  for 
ammonia  in  water,  the  contributions  of  the  different  elemen¬ 
tal  processes  in  Eq.  (17)  can  be  separated  and  inferences 
made  about  the  nature  of  the  species  formed  at  the  interface 
between  water  and  gaseous  ammonia. 

2.5.2.  Experimental  determination  of  the  mass 
accommodation  coefficient  a 

The  first  determinations  of  the  mass  accommodation  co¬ 
efficient  for  ammonia  were  made  indirectly  by  measuring 
the  evaporation  rates  of  ammonium  nitrate  and  of  chloride 
from  their  aqueous  solutions  [18-20].  In  the  mathematical 
treatment  of  the  experimental  data,  the  mass  accommodation 
coefficient  is  treated  as  an  adjustable  parameter,  selected  to 
fit  the  data.  In  these  experiments,  ammonia  evaporates  from 
the  solution  together  with  HCl  or  NO3H  and  possibly  with 
some  undissociated  salt.  The  values  of  a  reported  varied  be¬ 
tween  0.02  and  1. 

More  recent  determinations  of  a  by  Ponche  et  al.  [21], 
Bongartz  et  al.  [22],  and  Shi  et  al.  [9],  used  the  same 
common  approach  of  measuring  the  absorption  rate  of  am¬ 
monia  diluted  in  a  carrier  gas  into  a  liquid  jet  of  mono 
dispersed  droplets.  Results  of  these  three  investigations  are 
in  reasonable  agreement  when  differences  in  temperature 
and  measurement  uncertainties  are  taken  into  account,  as 
illustrated  in  Fig.  1.  The  mass  accommodation  coefficient 
increases  as  the  temperature  decreases,  from  0.08  at  290  K 
to  0.35  at  260  K. 

2.5.5.  Mechanism  of  ammonia  dissolution  in  water 

To  account  for  discrepancies  between  the  Gibbs  free 
energy  calculated  from  the  temperature  dependence  of  the 
accommodation  coefficients  of  a  series  of  alcohols  and  the 
Gibbs  free  energy  of  solvation  calculated  from  the  corre¬ 
sponding  thermodynamic  changes  of  enthalpy  and  entropy, 
Davidovits  and  co-workers  [23,24]  have  proposed  a  new 
model  for  the  uptake  of  gas  molecules  by  liquids.  The  model 
uses  the  concept  that  the  interface  is  a  narrow  region  of  a 
dense  gas-like  state  within  which  nucleation  is  constantly  oc¬ 
curring.  Incoming  gas  molecules  participate  in  the  formation 
of  molecular  clusters  between  solvent  and  gas  molecules. 
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Fig.  1.  TemperaUire  dependence  of  the  mass  accommodation  coefficient 
of  ammonia  in  water  [9]. 


When  the  cluster  reaches  a  critical  size,  it  is  incorporated 
in  the  bulk  of  the  liquid,  thereby  completing  the  solvation 
process  and  preventing  desorption  of  the  now  solubilized 
gas  molecules.  The  change  in  Gibbs  free  energy  calculated 
from  the  temperature  dependence  of  the  accommodation 
coefficient  is  interpreted  to  be  the  activation  energy  for  the 
formation  of  the  critical  size  cluster.  This  model  provides 
a  rational  explanation  for  a  surprising  series  of  measured 
changes  of  enthalpies  and  entropies  for  the  uptake  in  water 
of  vapor  of  alcohols  with  increasing  hydrophobicity. 

Donaldson  [25]  used  temperature-  and  time-dependent 
surface  tension  measurements  of  aqueous  ammonia  so¬ 
lutions  to  determine  the  interfacial  binding  energies  and 
evaporation  rates  of  ammonia.  The  values  for  enthalpy  and 
entropy  derived  from  the  equilibrium-tension  measurements 
are  consistent  with  almost  fully  solvated  ammonia  in  its 
surface-bound  state.  Ab  initio  calculations  suggest  that  am¬ 
monia  is  bound  to  a  small  number  (two  or  three)  of  water 
molecules  at  the  surface.  This  complex  species  would  rep¬ 
resent  the  “critical  cluster”  postulated  in  Davidovits  model 
that  would  be  easily  transferred  to  the  bulk  solution.  Spec¬ 
troscopic  analysis  by  Simonelli  and  Shultz  [26]  indicates 
that  ammonia  is  hydrogen  bonded  to  water  at  the  interface, 
but  that  ammonia  is  oriented  so  that  the  nitrogen  electron 
pair  is  bonded  to  the  water’s  OH  group  exposed  at  the 
liquid  surface. 

In  a  subsequent  publication,  McDuffie  [27]  analyzed 
surface  tension  and  vapor-liquid  equilibrium  data  over  the 
complete  range  of  ammonia-water  molecular  ratios.  Below 
a  molar  fraction  of  0.5  for  ammonia,  his  results  are  con¬ 
sistent  with  Donaldson’s  results  and  with  the  postulated 
ammonia-water  cluster.  Above  0.5  molar  fraction  of  am¬ 
monia,  an  abrupt  change  of  behavior  is  observed  that  would 
be  best  explained  by  the  formation  of  a  monolayer  of  NH3 
on  the  aqueous  solution. 


2  A.  Interactions  of  ammonia  with  atmospheric  water 


2A.L  In-cloud  scavenging  of  ammonia 
Shimshock  and  De  Pena  [28,29]  have  analyzed  the  mass 
transfer  of  gaseous  ammonia  to  a  raindrop  using  Slinn’s  [30] 
expression  for  the  flux  of  gas  to  a  drop  of  radius  R: 

^(W[Cgas])  =  4;r/?2A:o(^Pgas-^^^  (18) 


where  Cgas  is  the  concentration  of  the  dissolved  gas  in  the 
droplet,  pgas  the  partial  pressure  of  the  gas,  the  overall 
transfer  velocity  of  the  gas  in  both  gas  and  liquid  phase,  and 
a*  is  the  overall  solubility  parameter. 

In  the  case  of  anunonia,  they  assumed  an  accommoda¬ 
tion  coefficient  of  1  and  by  considering  the  dissociation  of 
ammonia  in  water  derived  the  following  expression  for  a*: 

a*  =  ^:H^l  +  [H30+]^^  (19) 


Because  of  the  high  solubility  of  ammonia  in  water,  they 
assumed  that  the  mass  transfer  in  the  liquid  will  be  much 
faster  than  the  mass  transfer  in  the  gas  phase  and  used  the 
following  approximate  equation  for  ATq: 


(20) 


where  D  is  the  diffusion  coefficient,  R  the  radius  of  the  drop, 
Sh  the  Sherwood  number  (=1  +  0.4(7?e)^/^(5c)^/^),  Re  the 
Reynolds  number,  Sc  the  Schmidt  number  (z=rjfD),  and  rj 
the  kinematic  viscosity  of  air. 

Integrating  Eq.  (18)  and  solving  to  obtain  the  time,  t, 
required  for  Cnhs  to  increase  from  an  initial  concentration 
[CNHali  to  a  final  concentration  [CNHalf  gives  Eq.  (21): 


\3Ko)  \Pnh,-  [CNHalf/aV 


(21) 


Using  Eq.  (21)  and  the  approximate  values  of  Kq  and  a* 
computed  from  Eqs.  (19)  and  (20),  the  time  to  reach  equilib¬ 
rium  concentration  of  ammonia  was  computed  for  different 
droplet  sizes. 

The  results  indicate  that  for  the  very  small  droplets,  R  = 
10-100  fxm,  characteristic  of  cloud  droplet  size,  the  equilib¬ 
rium  is  reached  in  a  matter  of  minutes,  shorter  than  the  length 
of  time  that  a  typical  droplet  exists  in  the  cloud.  Therefore, 
the  ammonia  in  cloud  droplets  will  be  in  equilibrium  with 
the  concentration  of  ammonia  in  the  interstitial  air  in  the 
cloud.  It  should  be  noted  that  this  conclusion  is  consistent 
with  the  conclusions  reached  through  alternative  approaches 
for  analysis.  For  example,  Ponche  et  al.  ([21]  and  references 
therein)  used  a  framework  of  analysis  for  mass  transport  lim¬ 
its  developed  by  Freiberg  and  Schwartz  [31]  and  concluded 
that  for  ammonia,  the  mass  transfer  process  will  be  limited 
by  liquid  phase  diffusion  in  the  droplet  and  that  the  time  to 
reach  equilibrium  for  a  100  jxm  droplet  will  be  0.32  s.  Al¬ 
though  the  estimates  of  time  are  very  different,  they  both 
indicate  that  for  the  small  diameter  droplets  found  in  clouds. 
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the  ammonia  concentration  in  the  droplet  will  be  in  equilib¬ 
rium  with  the  ammonia  partial  pressure  in  the  interstitial  air 
in  the  cloud. 

If  the  liquid  water  content  of  the  cloud  is  known,  the 
concentration  of  ammonia  in  the  droplet  can  be  calculated 
using  the  effective  Henry  law  constant  (Eq.  (19))  using  the 
algorithm  given  by  Shimshock  and  de  Pena  [28]. 

2.4,2,  BeloW‘Cloud  scavenging  of  ammonia 
For  the  larger  drop  size  characteristic  of  rain  showers, 
the  models  discussed  in  the  preceding  section  indicate  that 
the  concentration  of  ammonia  in  the  drop  will  not  have 
reached  equilibrium  with  the  ammonia  in  the  surrounding 
atmosphere  by  the  time  the  drop  reaches  the  ground.  Asman 
[32]  developed  the  most  rigorous  model  of  the  below-cloud 
scavenging  process.  The  model  calculates  the  concentration 
of  ammonia  in  the  raindrop  as  a  result  of  two  phenomena: 

(1)  Uptake  of  gaseous  ammonia  at  the  drop  surface  and 
subsequent  dissociation  and  liquid  diffusion. 

(2)  Change  in  drop  radius  due  to  water  evaporation. 

The  model  calculates  the  height  of  the  cloud  base  from 
dew  point  calculation  as  a  function  of  relative  humidity  and 
temperature  at  ground  level,  giving  Eq.  (22): 

Zbase  =  ci+^iln(RH(0))  (22) 

where  Cl  =  -9620+62. 10Ta(0),cfi  =  -5254-1 3. 43Ta(0), 
Ta(0)  is  the  temperature  at  ground  level  in  K,  RH(0)  the 
relative  humidity  at  ground  level  in  %,  and  Zbase  is  the  height 
of  cloud  base  in  meters. 

Asman  selected  Best’s  distribution  function  for  raindrop 
size  distribution  over  other  published  distributions  for  the 
following  reasons: 

•  Best’s  distribution  tends  to  give  the  smallest  scavenging 
coefficient.  Since  the  model  assumes  no  vertical  speed  of 
the  air  through  which  the  drops  fall,  the  calculated  scav¬ 
enging  coefficient  will  be  larger  than  the  actual  scaveng¬ 
ing  coefficient.  Choosing  the  distribution  that  gives  the 
smallest  estimate  of  the  scavenging  coefficient  therefore 
yields  a  value  closer  to  the  actual  value  under  real  atmo¬ 
spheric  conditions. 

•  The  Best  distribution  is  representative  of  observed  distri¬ 
bution  over  a  larger  range  of  rainfall  intensity  than  any 
other  model. 

Using  these  assumptions,  Asman  derived  Eq.  (23)  to  cal¬ 
culate  the  below-cloud  scavenging  coefficient  Xb: 

Xb  =  (23) 

where  I  is  the  rainfall  rate  at  ground  level,  mmh“*;  a:  aa  + 
bbDg  (ammonia  diffusivity  in  air),  aa:  oq  +  ai  RH(0);  bb: 
ho  +  b\  RH(0);  bav:  bavo  +  bavi  RH(0) 

With  the  following  values: 

•  ao  =  (4.476  x  IQ-^)  -  (1.347  x  lO-’)Ta(O); 

•  ai  =  (-3.004  X  10"’)  -  (1.498  x  10-®)Ta(0); 


Table  6 

Below-cloud  scavenging  coefficient  for  ammonia  (T  =  283  K;  RH  =  85%, 
^base  —  292  m) 


Rainfall  intensity  (mmh 

Scavenging  coefficient  (Asman)  (s  ^) 

1 

0.58  X  lO-* 

5 

1.60  X  10-* 

15 

3.1  X  10-^ 

25 

4.2  X  10-* 

•  feo  =  8.717  -  (2.787  x  10-2)Ta(0); 

•  bi=  (-5.074  X  10-2)  _  (2.894  x  10-4)Ta(0); 

•  bavo  =  (9.016  x  lO'^)  -  (2.315  x  10-2)Ta(0); 

•  bavi  =  (4.458  x  lO'^)  -  (2.115  x  10-5)Ta(0). 

The  amount  of  ammonia  scavenged  by  the  rainfall  can 
then  be  calculated  using  the  concentration  of  ammonia  in  the 
air  below  the  cloud,  volume  of  atmosphere  affected  by  the 
rain,  and  period  of  the  rainstorm,  knowing  the  temperature 
and  relative  humidity  at  ground  level. 

2,4,3,  Validation  and  application  of  the  models 

Comparison  of  predicted  scavenging  coefficients  with 
scavenging  coefficients  calculated  from  measured  concen¬ 
trations  of  ammonia  in  rainwater  has  been  less  than  satisfy¬ 
ing.  The  most  extensive  studies,  Sperber  and  Hameed  [33] 
and  Shimshock  and  de  Pena  [29],  fail  to  find  a  link  between 
measured  scavenging  coefficients  and  the  precipitation  in¬ 
tensity.  Both  groups  observed  a  wide  range  of  scavenging 
coefficients,  from  less  than  10“^  to  4.0  x  10""^  s~^  in 
Sperber’ s  study  and  from  6.6  x  10“^  to  8.5  x  10“^  s“^ 
in  Shimshock’ s  study.  On  the  positive  side,  the  median 
experimental  value  found  by  Shimshock  and  de  Pena, 
1.3  X  s"^  is  close  to  the  mean  value  found  by  Sperber 
and  Hameed,  1.2  x  10“^  s”^  (standard  deviation  0.2  x  10“"^). 
Both  are  in  the  same  range  as  values  predicted  by  Asman*  s 
model  (Table  6). 

Discrepancies  between  experimental  and  predicted  scav¬ 
enging  coefficients  can  be  credited  to  the  difficulty  of  achiev¬ 
ing  accurate  experimental  determinations  due  to: 

•  Difficulty  measuring  precisely  the  low  ammonia  concen¬ 
trations  normally  present  in  the  atmosphere  (a  few  ppbv). 

•  No  reliable  way  of  differentiating  between  ammonium 
ions  originating  from  ammonia  and  from  scavenged  am¬ 
monium  salts  in  collected  rainwater. 

•  Difficulty  measuring  actual  concentrations  of  ammonium 
salt-containing  particulates  in  the  atmosphere. 

•  Difficulty  separating  contributions  of  in-cloud  and 
below-cloud  scavenging  to  the  total  concentration  of 
anunonium  species  in  collected  rainwater. 

The  combination  of  these  factors  made  experimental  de¬ 
terminations  difficult  to  interpret  and  to  generalize.  This 
point  is  further  illustrated  by  the  wide  range  of  conclusions 
reached  by  the  researchers  estimating  the  proportion  of  the 
total  ammonia  deposition  contributed  by  below-cloud  scav- 
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enging:  Shimshock  and  de  Pena  [29]:  10%;  Meszaros  and 
Szentimrei  [34]:  26-74%;  and  Aneja  et  al.  [35]:  15%. 

In  summary,  the  Asman  model  appears  to  give  reasonable 
values  for  below-cloud  scavenging  and  can  be  used  to  esti¬ 
mate  scavenging  of  ammonia  by  raindrops  at  a  given  ground 
level  temperature  and  relative  humidity. 

2.5.  Interactions  with  the  surface  of  bodies  of  water:  *'dry 
deposition**  on  riven  lake,  and  sea  surfaces 

2,5.1.  Framework  for  analysis 

“Dry  deposition”  refers  to  the  process  by  which  atmo¬ 
spheric  pollutants  are  removed  from  the  atmosphere  with¬ 
out  significant  interaction  with  atmospheric  moisture  such  as 
wet  deposition:  removal  by  in-cloud  and  below-cloud  scav¬ 
enging,  fog  droplets,  or  dew.  For  water-soluble  gases  such  as 
ammonia,  dry  deposition  over  the  surface  of  natural  bodies 
of  water  (rivers,  lakes,  seas)  represents  an  important  mech¬ 
anism  for  removal  from  the  atmosphere. 

The  flux  of  removed  ammonia  over  a  unit  surface  area  of 
water  is  given  by  Eq,  (24): 

F  =  A[CNH3]gVdt  (24) 

where  [CisrHslg  is  the  concentration  of  ammonia  in  the  air  in 
jxgm”^;  for  NH3,  1  =  58.8nmol/w"^  =  1.32 ppbv 

(nmol  =  nanomole  =  10”^  mol),  Vd  the  dry  deposition 
velocity  in  ms~^,t  is  duration  of  the  event  in  s,  F  the  flux 
of  deposited  ammonia  in  kg  m“^,  and  A  the  unit  conversion 
constant  (=10’”^). 

The  dry  deposition  velocity  is  really  an  expression  of  the 
mass  transfer  coefficient  of  ammonia  from  the  air  into  the 
water  phase.  Two  general  approaches  have  been  used  by 
Shahin  et  al.  [36]  to  model  this  mass  transfer  process. 

2.5. 7. 7.  Application  of  the  film  theory.  The  film  theory 
posits  that  the  overall  mass  transfer  resistance  is  imposed 
over  a  thin  film  at  the  interface  or  at  the  boundary  layer 
through  which  diffusion  is  the  driving  force.  It  follows  that 
the  mass  transfer  coefficient  is  directly  proportional  to  the 
diffusion  coefficient  and  to  the  thickness  of  the  boundary 
layer.  In  the  case  of  the  surface  of  a  natural  body  of  water, 
the  thickness  of  the  boundary  layer  is  controlled,  to  a  large 
extent,  by  the  wind  speed.  Several  empirical  relationships 
between  the  mass  transfer  coefficient  and  the  wind  speed 
have  been  developed  for  water-soluble  gases  such  as  am¬ 
monia.  Based  on  measurements  of  water  evaporation  rate, 
Thibodeaux  [37]  derived  more  elementary  but  relatively 
complex  expressions  for  the  mass  transfer  coefficients  for 
laminar  and  for  turbulent  flows.  A  critical  review  of  this 
approach  is  presented  in  Shahin  et  al.’s  paper  [36]. 

2.5. 7.2.  Application  of  the  resistance  model.  The  resis¬ 
tance  model  is  the  most  widely  used  model  for  atmospheric 
deposition  processes  [36,38-40].  By  analogy  with  electrical 
resistances,  the  deposition  velocity,  Vd»  is  inversely  propor¬ 


tional  to  the  sum  of  three  resistances  in  series: 
Vd=  ^ 


Ra  -|-  Rb  -|-  Rch 


(25) 


where  Ra  is  theaerodynamic  resistance,  i.e.,  resistance  of 
transport  through  the  turbulent  surface  layer,  Rb  the  resis¬ 
tance  for  diffusion  through  quasi-laminar  air  layer  in  con¬ 
tact  with  the  water  surface,  and  Rch  the  resistance  due  to 
physicochemical  phenomena  occurring  at  the  interface. 

The  aerodynamic  resistance,  Ra,  depends  on  wind  speed, 
atmospheric  stability,  and  surface  roughness.  Typically,  Ra 
is  calculated  from  the  wind  velocity,  at  the  reference 
height  at  which  the  concentration  of  ammonia  is  measured 
(typically  10  m  above  the  surface),  and  from  the  friction 
velocity,  m*: 


Ra: 


(W*)2 


(26) 


Several  semi-empirical  expressions  have  been  developed 
to  estimate  the  friction  velocity  from  the  roughness  length  of 
the  surface,  if  known,  or  from  the  standard  deviation  of  the 
wind  direction.  Shahin  et  al.  [36]  used  this  previous  approach 
that  appears  best  suited  for  incorporation  in  a  dispersion 
model.  The  aerodynamic  resistance  is  then  given  by: 


Ra=: 


- * ,  for  neutral  and  stable  atmospheres 

UIOOq 

9 

- * ,  for  unstable  atmospheres 

uioa^ 


(27) 


where  wio  is  the  wind  velocity  10  m  above  the  surface  and 
ae  is  the  standard  deviation  of  the  wind  direction  in  radians. 

The  boundary  layer  resistance,  Rb,  is  a  function 
of  the  molecular  diffusion  coefficient  for  ammonia,  D 
(2.09  X  10”^  m^  s"”^  at  283K),  of  the  kinematic  viscosity  of 
air  (v  =  1.47  x  lO”^  m^  s“^  at  283  K),  and  of  the  friction 
viscosity,  m*,  and  can  be  estimated  by  using  one  of  several 
published  semi-empirical  relationships,  such  as  Eq.  (28): 

Rb=_(^)'  (28) 

ku*  \DPr) 


where  k  is  the  von  Karman’s  constant  (0.4)  and  Pr  is  the 
Prandt  number  of  air  (0.72). 

In  general  Rch  is  measured  experimentally  from  labora¬ 
tory  experiments  using  static  chambers. 


2.5.2.  Determination  of  the  deposition  velocity 

Three  major  approaches  have  been  used  to  measure  the 
deposition  velocity  of  ammonia  over  the  surface  of  natural 
bodies  of  water: 

(1)  Laboratory  determination  of  Rch  and  calculations  of  Ra 
and  Rb  as  outlined  in  the  preceding  section  [38]. 

(2)  Field  measurements  of  ammonium  compound  deposi¬ 
tion  over  extended  period  of  time  [41,42]. 

(3)  Field  measurements  using  Water  Surface  Sampler 
(WSS)  and  Knife-edge  Surrogate  Surface  (KSS)  tech¬ 
niques  [36]. 
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Table  7 


Deposition  velocity  of  ammonia  over  seawater 


Location 

Vd  (mms  *) 

Reference 

Atlantic  Basin 

8-20 

Quinn  et  al.  [44] 

Australia-Southern  Ocean 

8.3 

Griffiths  et  al.  [45] 

NE  Pacific 

8.3 

Quinn  et  al.  [46] 

North  and  Baltic  Seas 

8 

Barrett  [47] 

North  Sea 

7.6 

Asman  et  al.  [48] 

Tampa  Bay 

7 

Poor  et  al.  [42] 

North  Sea 

2-15 

Lee  et  al.  [38] 

Larsen  et  al.  [43]  have  compiled  all  the  results  published 
so  far  for  the  deposition  velocity  of  ammonia  over  seawater 
(Table  7). 

Although  the  results  in  Table  7  appear  to  be  relatively 
consistent,  the  values  reported  are  averages  of  measure¬ 
ments  with  large  uncertainties  associated  with  them.  When 
reported,  the  range  of  values  quoted  reflects  the  uncertainty 
of  the  measurements  due  principally  to  the  difficulty  of  mea¬ 
suring  very  small  concentrations  of  ammonia  in  the  air  and 
to  the  lack  of  easy  differentiation  between  ammonia  and 
ammonium  species  from  particles  when  dissolved  in  water. 
These  values  should  be  regarded  as  order  of  magnitude  only. 

Shahin  et  al.  [36]  conducted  extensive  measurements 
and  analysis  of  dry  deposition  of  ammonia,  nitric  acid,  and 
sulfur  dioxide  over  water  surfaces  in  Chicago.  They  used 
direct  measurements  of  dry  deposition  using  surrogate  tech¬ 
niques.  An  aerodynamically  designed  Water  Surface  Sam¬ 
pler  collects  the  total  flux  of  the  species  considered  while 
a  Knife-edge  Surrogate  Surface  sampler  captures  only  the 
particulates  from  an  identical  air  stream  sample.  The  dif¬ 
ference  between  these  two  measurements  gives  an  accurate 
measurement  of  the  true  dry  deposition  of  a  water-soluble 
gas  such  as  ammonia.  Their  combined  data  for  ammonia, 
nitric  acid,  and  sulfur  dioxide  fits  the  following  model: 

Vd  =  D®*5{(0.98  ±  O.Dmio  +  (1.26  =b  0.3))  (29) 

where  D  is  the  gas-phase  diffusion  coefficient  of  the 
water-soluble  gas,  here  ammonia,  and  wio  is  the  wind  speed 
10  m  above  the  water  surface;  the  ±values  are  the  95% 
confidence  intervals  for  the  corresponding  constants.  Tem¬ 
perature  was  found  to  have  little  effect  on  Comparison 
with  the  predictions  of  Thibodeaux’s  film  theory  model 
found  that  Thibodeaux’s  model  underestimates  Vd  while 
the  empirical  model  based  on  the  resistance  model  tends  to 
slightly  overestimate  Vj  (see  Table  2  in  Shahin’ s  paper). 

The  values  predicted  by  Eq.  (29)  (Table  8)  are  in  gen¬ 
eral  higher  than  the  average  values  reported  for  seawater  in 
Table  8.  A  possible  explanation  would  be  that  Shahin  et  al. 
[36]  used  fresh  water  in  their  WSS  apparatus  where  ammo¬ 
nia  should  be  expected  to  be  more  soluble  than  in  seawater 
(salting-out  effect). 

2.5.J.  Recommendations 

The  dry  deposition  of  ammonia  over  natural  bodies  of 
water  is  likely  to  be  a  significant  mechanism  for  the  removal 


Table  8 

Deposition  velocity  of  ammonia  on  fresh  water  (Eq.  (28)) 


Wind  speed,  wio  (ms  ^) 

Deposition  velocity,  Va  (mms  ^) 

1 

11.4 

2 

16.4 

3 

21.4 

4 

26.4 

5 

31.4 

6 

36.4 

of  ammonia  from  the  atmosphere.  The  data  presented  in 
the  preceding  sections  have  been  determined  for  the  very 
low  concentrations  of  ammonia  normally  present  in  the 
atmosphere,  i.e.,  0.1“5ppbv.  In  the  case  of  accidental  re¬ 
lease  of  liquid  ammonia  from  a  tank  car,  the  concentration 
of  ammonia  in  the  air  will  be  several  orders  of  magnitude 
greater  than  normally  present  in  the  atmosphere.  The  driv¬ 
ing  forces  and  rate  limiting  steps  for  the  dry  deposition  of 
ammonia  over  water  surfaces  will  be  drastically  different 
from  those  prevailing  under  normal  atmospheric  conditions. 
Therefore,  it  is  recommended  that  laboratory  experiments 
be  conducted  to  estimate  the  deposition  velocities  of  high 
gas-phase  concentrations  of  ammonia  over  both  fresh  water 
and  simulated  seawater. 


3.  Dry  deposition  of  ammonia  over  the  vegetation 

3,L  Mechanism  of  uptake  of  ammonia  by  plants 

Dry  deposition  is  a  mass  transfer  process  whereby  ammo¬ 
nia  is  first  transported  to  the  surface  by  turbulent  and  molec¬ 
ular  diffusion  and  then  removed  by  adsorption  or  absorption 
at  the  surface.  In  the  case  of  vegetation,  die  situation  is  fur¬ 
ther  complicated  by  the  role  of  ammonia  in  plant  physiology. 
Depending  on  the  concentration  of  ammonia  prevailing  in 
the  atmosphere,  vegetation  can  be  either  a  source  or  a  sink 
for  ammonia.  This  two-way  exchange  of  ammonia  between 
vegetation  and  atmosphere  has  become  better  understood  in 
the  past  few  years,  reconciling  numerous  contradictory  re¬ 
sults  reported  in  earlier  publications. 

In  1980,  Farquhar  et  al.  [49]  suggested  growing  vege¬ 
tation  would  emit  ammonia  when  exposed  to  air  with  a 
concentration  of  NH3  below  some  “compensation  point” 
related  to  the  partial  pressure  of  NH3  in  the  substomatal 
cavities.  Subsequent  studies  confirmed  the  existence  of  this 
compensation  point  and  found  its  value  was  in  the  range 
of  0.5-4  fxg  m”^  (0.66-5.28  ppbv)  with  a  median  value  be¬ 
tween  1.5  and  2|xgm”^  (2-2.64  ppbv)  [50-56].  Distinction 
can  be  made  [57]  between  stomatal  and  canopy  compensa¬ 
tion  points.  The  stomatal  compensation  point  refers  solely 
to  the  adsorption  through  the  stomata,  the  minute  openings 
in  the  epidermis  of  a  plant  through  which  gaseous  exchange 
with  the  atmosphere  takes  place.  Canopy  compensation  point 
includes  the  potential  exchange  of  ammonia  between  the 
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whole  plant  surface  and  the  atmosphere,  considering  that 
particularly  at  elevated  humidity  (RH  >  70%)  dew  on  the 
leaf  surface  can  adsorb  ammonia  reversibly  or  that  the  waxy 
surface  of  the  leaves  (cuticle)  can  retain  ammonium  salt  par¬ 
ticles  (nitrate,  chloride)  that  are  in  chemical  equilibria  with 
the  vapor  of  their  constituents.  Van  Hove  et  al.  [51]  did 
an  extensive  study  of  the  stomatal  compensation  point  for 
rye  grass  and  concluded  this  value  varied  between  0.5  and 
4jxgm“^  and  was  influenced  mainly  be  temperature.  The 
typical  concentration  of  ammonia  varying  from  0.75  p.g  m“^ 
in  pristine  troposphere  to  as  much  as  19  fxgm”^  in  polluted 
air  ([40]  p.  37),  the  stomatal  compensation  point  is  generally 
exceeded  and  the  vegetation  acts  usually  as  a  sink  for  atmo¬ 
spheric  ammonia.  This  will  definitely  be  the  case  during  a 
release  incident. 


Table  9 

Land-dependent  values  for  canopy  resistance  to  ammonia  mass  transfer 


Land  type 

Canopy  resistance  (smm 

Arable  soil 

1 

Forest 

0.02 

Grass 

0.6 

Moorland 

0,02 

Urban 

0.24 

laboratory  and  field  dry  deposition  measurements  for  reac¬ 
tive  nitrogen  compounds,  including  ammonia. 

3.3 .L  Leaf -level  measurements 

Three  common  approaches  are  used  in  determining  the 
mass  transfer  coefficient  of  ammonia  to  exposed  vegetation: 


3.2.  Analysis  and  modeling  of  the  dry  deposition  process 


A  resistance  model  similar  to  the  one  described  in  the  pre¬ 
ceding  section  for  water  surfaces  is  most  commonly  used  to 
analyze  the  dry  deposition  process.  The  deposition  velocity, 
Vd»  is  inversely  proportional  to  the  sum  of  three  resistances 
in  series  equation: 


Vd  = 


1 

Ra  +  Rb  +  Rc 


(25) 


where  Ra  is  the  aerodynamic  resistance,  i.e.,  resistance  of 
transport  through  the  turbulent  surface  layer,  Rb  the  resis¬ 
tance  for  diffusion  through  quasi-laminar  air  layer  in  con¬ 
tact  with  the  soil  or  vegetation  surface,  and  Rc  the  canopy 
resistance. 

Aerodynamic  and  laminar  resistances  have  the  same 
meaning  and  mathematical  description  as  given  in  Section 
2.5.  The  canopy  resistance  represents  a  composite  of  the  var¬ 
ious  physico-chemical  phenomena  involved  in  the  interac¬ 
tion  of  ammonia  with  the  surface.  This  includes  the  ammonia 
uptake  through  the  stomata  of  leaves,  adsorption  of  ammonia 
by  moisture  of  the  impacted  surfaces  (leaves,  bark),  as  well 
as  the  absorption  and  chemical  transformation  of  ammonia 
with  soil  components  [58,59].  Because  of  the  roughness 
of  the  vegetation,  canopy  resistance  generally  dominates 
the  process,  particularly  during  the  daytime.  It  is  seldom 
possible  to  obtain  separate  values  for  the  three  component 
resistances  of  the  deposition  velocity  and  experimentally  de¬ 
termined  values  for  are  generally  used  (see  next  section). 
However,  in  some  large  scale  modeling  efforts,  when  a  great 
variety  of  terrains  had  to  be  considered,  generally  accepted 
average  values  for  Rc  (Table  9)  have  been  used  [60]. 


3.3.  Experimental  determination  of  dry  deposition  velocity 

Measurement  of  dry  deposition  velocities  are  usually  con¬ 
ducted  at  two  levels,  individual  plant  leaves  in  laboratory 
settings  and  forest  canopy  in  the  field  level.  Hanson  and 
Lindberg  [61]  published  a  comprehensive  review  of  both 


(1)  Mass  balance:  vegetation  (leaves,  branches,  whole 
plant)  is  exposed  to  a  set  concentration  of  ammonia  in 
air  in  an  environmentally  controlled  chamber  and  am¬ 
monia  concentration  is  measured  at  the  chamber  inlet 
and  outlet  to  determine  ammonia  uptake  of  the  plant. 

(2)  Isotopic  labeling:  plant  is  exposed  to  ^^NHa  and  the 
level  of  radioactive  label  incorporated  into  the  plant  tis¬ 
sue  is  measured. 

(3)  Foliar  extraction:  concentration  of  nitrogen  species  is 
determined  in  an  environmentally  controlled  chamber 
in  plant  tissue  before  and  after  ammonia  exposure. 

The  mass  balance  technique  is  the  most  widely  used  of 
the  three  approaches. 

3.3.2.  Whole  canopy  measurements 

The  whole  canopy  measurement  approach  uses  microme- 
teorological  methods  based  on  an  assumption  that  the  verti¬ 
cal  flux  through  a  horizontal  plane,  measured  at  a  reference 
level  above  the  surface,  provides  the  flux  through  the  at¬ 
mosphere/vegetation  interface.  In  general,  two  common  ap¬ 
proaches  are  normally  used: 

(1)  The  eddy  correlation  technique:  The  eddy  correlation 
technique  measures  vertical  turbulent  flux  directly  from 
calculations  of  the  mean  covariance  between  the  verti¬ 
cal  component  of  the  wind  and  between  the  pollutant 
concentration.  Since  fast  and  accurate  ammonia  concen¬ 
tration  measurement  techniques  are  not  available,  heat 
fluxes  and  local  ammonia  concentration  measurements 
are  used  to  approximate  the  actual  turbulent  flow  mea¬ 
surements. 

(2)  The  flux  gradient  or  profile**  technique:  Using  the  flux 
gradient  or  “profile”  technique,  the  flux  is  derived  from 
measurements  of  the  vertical  concentration  profile  of 
ammonia  and  eddy  exchange  coefficients. 

Because  of  the  various  difficulties  involved  in  measuring 
in  real  time,  extremely  low  concentrations  of  ammonia  in  the 
atmosphere,  the  whole  canopy  measurements  for  ammonia 
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Table  10 

Mass  transfer  coefficient  of  NH3  to  leaf  surfaces  (chamber  method,  day¬ 
time  conditions) 


Plant  species 

NH3  cone. 

(ixgm'3) 

Conductance, 
1/Rc  (mms”^) 

Reference 

Fescue 

341 

15.3 

Rogers  and  Aneja  [63] 

Soybean 

37 

4.4 

Hutchinson  et  al.  [64] 

170 

11 

Rogers  and  Aneja  [63] 

Cotton 

63 

2.2 

Hutchinson  et  al.  [64] 

331 

6.7 

Rogers  and  Aneja  [63] 

Sunflower 

45 

4.3 

Hutchinson  et  al.  [64] 

Orchard  grass 

283 

10.3 

Rogers  and  Aneja  [63] 

Snap  bean 

8 

6-32 

Farquhar  et  al.  [49] 

140 

13 

Rogers  and  Aneja  [63] 

144 

2-6 

van  Hove  et  al.  [50] 

502 

2-6 

van  Hove  et  al.  [50] 

Poplar 

143 

0.5-9 

van  Hove  et  al.  [50] 

Wheat 

277 

15.1 

Rogers  and  Aneja  [63] 

Com 

34 

6.5 

Hutchinson  et  al.  [64] 

320 

3.6 

Rogers  and  Aneja  [63] 

have  generally  been  accomplished  using  the  flux  gradient 
method. 

33,3.  Results  obtained 

The  reported  values  for  the  dry  deposition  velocity  of 
ammonia  are  listed  in  Table  9  for  whole  canopy  experi¬ 
ments  and  in  Table  10  for  leaf  surface  measurements.  Al¬ 
though  they  are  expressed  in  the  same  units,  the  values  re¬ 
ported  for  leaf  surface  and  for  whole  canopy  experiments 
are  not  directly  comparable  because  they  involved  differ¬ 
ent  receptor  areas  and  characterize  the  process  at  a  different 
level  of  substitution.  It  is  possible  to  convert  the  mass  trans¬ 
fer  coefficient  obtained  on  leaf  surface  measurements  to  a 
true  deposition  velocity  by  multiplying  it  by  an  appropriate 
leaf  area  index  (LAI)  that  is  a  representative  characteristic 
of  the  plant  species  considered.  The  theoretical  derivation 
of  this  LAI  is  complex  but  semi-empirical  values  for  vari¬ 
ous  plant  species  have  been  published  (e.g.,  Hanson  et  al. 
[62]). 

An  inspection  of  the  values  reported  in  Tables  10  and  1 1 
leads  to  the  following  conclusions: 


•  Deposition  velocities  ranging  from  about  4mms"^  to 
about  15mms“'^  characterize  annual  plants  of  modest 
height  and  foliage  crown  development. 

•  Deposition  velocities  ranging  from  about  25-50  mm  s**^ 
characterize  forests. 

•  There  is  a  large  variability  in  the  data  that  can  be  attributed 

o^’differences  of  experimental  techniques  and  data  analy¬ 
sis  methodology  among  investigators; 
o  difficulties  of  measuring  accurately  low  atmospheric 
concentrations  of  ammonia; 

o  influence  of  determining  variables  not  controlled  or  av¬ 
eraged  out  in  the  experimental  design  e.g.,  humidity, 
temperature,  sunlight. 

In  the  most  recent  investigations,  the  influence  of  these 
external  factors  has  been  quantified  and  found  to  have  a 
significant  impact  of  the  measured  deposition  velocities  (see 
next  section). 

3,3,4,  Factors  influencing  deposition  velocities 

Influence  of  humidity:  It  has  long  been  suspected  that 
condensed  atmospheric  moisture  deposited  on  the  surface  of 
leaves  constituted  a  potential  sink  for  atmospheric  ammonia 
[72-74].  Wyers  and  Erisman  [75]  have  attempted  to  quan¬ 
tify  the  impact  of  the  moisture  deposited  on  leaf  surfaces 
by  considering  the  Canopy  Water  Storage  (CWS)  level  of 
the  forest.  This  value  can  be  calculated  from  microwave 
measurement  of  the  water  content  of  the  canopy.  A  canopy 
saturated  with  water  (e.g.,  after  a  rainfall)  has  a  CWS  value 
exceeding  2  mm.  A  dry  canopy  has  a  CWS  value  of  less 
than  0.25  mm.  In  between,  the  canopy  is  drying.  When 
the  canopy  is  saturated  with  water  (CWS  >  2  mm)  the  re¬ 
sistance  to  ammonia  deposition  is  very  low  and  therefore 
the  deposition  velocities  are  at  their  maximum,  averag¬ 
ing  about  60mms“^  When  the  canopy  is  dry  (CWS  < 
0.25  mm),  the  deposition  velocity  decreases  exponentially 
with  the  decreasing  humidity  level.  Sutton  et  al.  have  de¬ 
rived  Eq.  (30)  from  both  field  and  laboratory  measure¬ 
ments: 

Vd  =  (50mms"^)  x  ^exp  - 


Table  11 

Deposition  velocity  for  whole  canopy 


Plant/terrain 

NH3  cone.  (|xgm~^) 

Vd  (mms^^) 

Method 

Reference 

Alpine  Tundra 

13  ±  6 

Gradient 

Rattray  and  Sievering  [65] 

Spruce  forest 

0.47 

26  (-125  to  201) 

Gradient 

Andersen  and  Hovmand  [66] 

Alpine  Timdra 

3.8 

4 

Estimate 

Langford  and  Fehsenfeld  [56] 

Douglas  Fir  Forest 

6.1 

22  (S,D.  50) 

Eddy  correlation 

Duyzer  et  al.  [67] 

Douglas  Fir  Forest 

0.1-25 

32  (median) 

Gradient 

Wyers  et  al.  [68] 

Douglas  Fir  Forest 

36  (24  h  average)  max  50 
midday  min  28  night 

Eddy  correlation 

Duyzer  et  al.  [69] 

Heathland 

16 

Eddy  correlation 

Duyzer  et  al.  [70] 

Soybean 

59 

6.1 

Chamber 

Aneja  et  al.  [71] 

Snap  bean 

76 

4 

Chamber 

Aneja  et  al.  [71] 

Fescue 

456 

11.6 

Chamber 

Aneja  et  al.  [71] 

Com 

190 

3.1 

Chamber 

Aneja  et  al.  [71] 
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Table  12 


Annual  median  deposition  velocities  during  nighttime  for  ammonia  on 
the  dry  canopy  of  coniferous  forest 


Year 

Annual  median  deposition  velocity  (mms 

RH  =42% 

RH  =  95% 

1993 

6 

100 

1994 

8.5 

32 

For  a  drying  canopy  (0.25  mm  <  CWS  <  2  mm)  the  average 
deposition  velocity  would  range  from  50  to  60mms“^ 

These  average  values,  extracted  from  Wyers  and  Erisman 
data,  are  only  applicable  to  well  developed  coniferous  for¬ 
est  stands.  For  terrains  with  a  lower  aspect  ratio  such  as 
pasture,  agricultural  fields,  wasteland  with  low  vegetation, 
the  deposition  velocities  are  probably  about  one-third  of  the 
corresponding  values  for  forests,  according  to  the  general 
trend  apparent  in  Tables  10  and  11. 

Daytime/nighttime:  Deposition  of  ammonia  on  vegeta¬ 
tion  during  the  night  is  influenced  by  both  physiological 
and  meteorological  characteristics  of  nighttime  compared 
to  daytime.  During  the  night,  the  stomata  of  plants  close 
thereby  stopping  ammonia  uptake  by  plant  tissues.  From 
a  meteorological  standpoint,  calmer  air,  lower  temperature 
and  therefore  lower  humidity  during  the  night  also  affect 
the  effectiveness  of  ammonia  deposition  when  night  and 
day  situations  are  compared.  For  coniferous  forests,  Duyzer 
et  al.  [69]  have  found  a  clear  maximum  for  the  deposition 
velocity  of  ammonia  at  midday  (50  mm  s"^ )  and  a  minimum 
of  about  28nuns“^  during  the  night.  Wyers  and  Erisman 
[75]  documented  a  similar  trend.  These  authors  also  quan¬ 
tified,  for  nighttime  conditions,  the  exponential  decrease  of 
deposition  velocity  with  decreasing  humidity  found  by  Sut¬ 
ton  for  daytime  conditions  (Eq.  (30)).  However,  their  values 
are  substantially  different  from  those  reported  by  Duyzer 
et  al.  and  the  annual  median  deposition  velocities  are  also 
variable  for  the  two  years  analyzed  (Table  12). 

Temperature  and  atmospheric  concentration  of  ammonia: 
In  field  experiments,  the  variability  of  the  ammonia  concen¬ 
tration  is  generally  relatively  small  and  its  effect  difficult  to 
isolate.  Similarly,  it  is  difficult  to  separate  the  effect  of  tem¬ 
perature  from  the  effect  of  humidity.  However,  under  the 
well-controlled  conditions  of  whole  canopy  chamber  exper¬ 
iments,  Aneja  et  al.  [71]  have  found  that  the  concentration 
of  ammonia  in  the  gas  phase  (from  90  to  1400  ppbv)  and 
the  temperature  (from  12  to  30  °C)  have  a  negligible  impact 
on  the  deposition  velocity.  However,  at  low  temperatures, 
the  decrease  in  physiological  activities  in  plants  leads  to  a 
closing  of  the  stomata  and  to  a  corresponding  increase  in 
resistance  to  ammonia  uptake  which  result  in  deposition  ve¬ 
locities  below  5  mms*"^  [73]. 

3.4.  Conclusions — recommendations 

The  dry  deposition  of  ammonia  over  vegetation  is  a  com¬ 
plex  phenomenon  that  is  still  poorly  understood  and  difficult 


to  determine  experimentally.  In  the  case  of  a  massive  re¬ 
lease  of  ammonia  in  the  atmosphere,  the  high  concentration 
of  ammonia  will  become  the  main  driving  force  for  the 
deposition  process.  Aerodynamic  and  gas-phase  diffusion 
resistance  will  become  negligible  and  the  physico-chemical 
phenomena  on  surfaces  will  control  deposition  velocity. 
Given  the  uncertainties  associated  with  the  results  in  the 
literature  and  the  difficulties  of  conducting  meaningful  ad¬ 
ditional  experiments,  we  recommend  the  adoption  of  the 
following  values  for  dry  deposition  velocities: 

•  For  forested  land: 

o  maximum  deposition  velocity  =  60mms”^  under  wet 
conditions  and  at  100%  humidity; 
o  apply  Sutton’s  exponential  decrease  (Eq,  (30))  to  calcu¬ 
late  the  deposition  velocity  for  lower  relative  humidity 
levels; 

o  for  nighttime  conditions,  the  maximum  deposition  ve¬ 
locity  (100%  humidity)  =  28mms"^ 

•  For  terrains  with  low  aspect  ratio  vegetation:  agricultural 
crops,  pastures,  tundra,  wasteland: 

o  maximum  deposition  velocity  =  20mms“^  under  wet 
conditions  and  at  100%  humidity; 
o  influence  of  humidity  level:  use  Sutton’s  exponential 
decrease  for  maximum  deposition  velocity  value, 
o  for  nighttime  conditions  assume  a  maximum  deposition 
velocity  value  of  10mms“^  at  100%  humidity. 


4.  Reactions  of  ammonia  with  atmospheric  sulfur  oxides 

4.1.  Reactions  with  sulfur  dioxide 

4.1.1.  Products  of  reaction 

•  Anhydrous  reactions  ofNH'i  and  SO2'  For  over  a  century 
and  a  half  chemists  have  known  that  these  two  gases  re¬ 
act  readily  but  which  products  are  formed  has  remained  a 
subject  of  controversy.  Becke-Goehring  [76]  and  Meyer 
et  al.  [77]  conducted  a  comprehensive  summary  of  early 
work.  A  more  coherent,  yet  complex  picture  emerges  from 
recent  work,  once  the  importance  of  the  molar  ratio  of 
the  two  gases  is  realized.  Vapor  pressure  measurements 
[78,79]  and  spectroscopic  analysis  [77]  lead  to  the  fol¬ 
lowing  conclusions: 

o  At  a  stoichiometric  ratio  of  1:1,  NH3  and  SO2  reacts 
according  to  Becke-Goehring  [76]  and  Meyer  et  al.  [77] 
according  to  Eq.  (31)  to  give  an  adduct  in  the  form  of 
a  yellow  solid,  in  equilibrium  with  the  vapor  pressure 
of  its  components: 

NH3  +  SO2  NH3SO2  (31) 

LH°  =  32.2  kcal  mol"  ‘ ;  AS°  =  84.8  cal  moP  ^  K"' 

o  At  a  stoichiometric  ratio  of  2:1,  NH3  and  SO2  reacts 
according  to  Eq.  (32)  to  give  an  adduct  in  the  form  of 
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a  white  solid,  in  equilibrium  with  the  vapor  pressure  of 
its  components: 

2NH3  +  SO2  (NH3)2S02  (32) 

A//°  =  62.2kcalmol“^;  AS°  =  174.8 cal mol“^  K“^ 

The  reactions  of  anhydrous  SO2  and  NH3  gases  lead¬ 
ing  to  the  formation  of  NH3*S02  and/or  (NH3)2*S02 
are  completely  reversible.  When  the  gas  vapor  pres¬ 
sure  is  sufficiently  low,  the  solids  dissociate  into  their 
original  components:  SO2  and  NH3.  The  equilibrium 
conditions,  i.e.,  the  pressure  at  which  the  solids  ap¬ 
pear  are  a  function  of  temperature.  For  the  formation 
of  NH3-S02,  the  product  of  the  partial  pressures  of  the 
reactants  at  equilibrium  varies  from  48Torr^  at  5°C 
(278. 15  K)  to  425  Torr^  at  45  °C  (3 1 8. 15  K).  For  the  for¬ 
mation  of  (NH3)2*S02  the  temperature  dependence  is 
steeper  and  the  corresponding  equilibrium  constant  in¬ 
creases  from  225  Torr^  at  5  °C  (278. 15  K)  to  1030  Torr^ 
at  15  (288.15  K)  [78]. 

o  In  the  presence  of  a  stoichiometric  excess  of  SO2, 
c/5-thionylimide,  HNSO,  is  initially  formed  at  low  tem¬ 
peratures,  followed  by  the  formation  of  red  heteroge¬ 
neous  (HNSO)jc  polymers  that  upon  warming  undergo 
auto-redox  decomposition  yielding  mixtures  of  sulfate, 
sulfamic  acid,  dithionate,  thiosulfate,  polythionate,  and 
elemental  sulfur.  For  stoichiometric  ratios  between  1:1 
and  2:1,  mixtures  of  the  two  adducts  are  obtained,  a 
situation  that,  together  with  the  unrecognized  influence 
of  traces  of  water  or  oxygen,  created  a  lot  of  confusion 
in  early  investigations. 

•  Influence  of  water  on  the  reaction  of  NH3  and  SO2: 
In  1975,  Heicklen  and  co-workers  [78,80]  elucidated 
the  critical  role  played  by  water  in  the  distribution  of 
products  obtained  from  the  reaction  of  NH3  with  SO2. 
If  water  is  present  in  excess  of  the  required  stoichio¬ 
metric  amount  (Eq.  (33)),  both  the  1:1  and  2:1  adducts 
react  quickly  to  form  stable  ammonium  sulfite.  When 
sub-stoichiometric  levels  of  water  are  present  ammonium 
pyrosulfite,  (NH4)2S205,  is  formed. 

2NH3(g)  +  S02(g)  +  H20(g)  (NH4)2S03(s)  (33) 

AH°  =  -53  kcal  mol"  ‘ ;  A5°  =  - 1 13  cal  moP  * 

4.7.2.  Oxidation  of  SO2  in  the  presence  of  ammonia 
In  1958,  Junge  and  Ryan  [81]  observed  that  the  rate  of 
oxidation  of  SO2  by  oxygen  was  faster  when  the  experiment 
was  conducted  using  aerosol  formed  from  sodium  chloride 
solution  spray  in  a  polyethylene  bag  than  when  the  reac¬ 
tion  was  conducted  in  bulk  aqueous  solutions.  The  same  ob¬ 
servation  was  made  when  ammonia  was  added  to  the  bag 
containing  a  water  aerosol  and  SO2.  Subsequent  investiga¬ 
tions  of  the  oxygen  oxidation  of  SO2  in  the  presence  of  wa¬ 
ter  aerosols  containing  sodium  chloride  or  sea  salts  [82,83] 
confirmed  these  observations.  Van  Den  Heuvel  and  Mason 
[84]  re-examined  the  oxidation  of  SO2  in  the  presence  of 


ammonia  and  concluded  that  the  mass  of  ammonium  sulfate 
formed  was  proportional  to  the  product  of  the  surface  area  of 
the  water  drops  and  the  time  of  exposure.  Although  the  po¬ 
tential  catalysis  of  the  reaction  by  trace  amounts  of  transition 
metal  ions,  particularly  manganese,  could  not  be  completely 
discounted  [85,86],  these  authors  concluded  that  the  acceler¬ 
ation  of  the  oxidation  reaction  was  attributable  to  the  special 
conditions  existing  at  the  surfac6  of  the  liquid  droplets. 

Brenner  and  co-workers  [87,88]  conducted  a  thorough, 
carefully  designed  investigation  of  the  oxidation  of  SO2  in 
the  presence  of  ammonia  under  simulated  cloud  conditions 
and  in  thin  water  films.  They  established  convincingly  that 
metal  catalysis  was  not  the  reason  for  the  increased  reaction 
rate.  Nevertheless,  when  anunonia  was  present,  80%  of  the 
SO2  was  oxidized  to  ammonium  sulfate  in  less  than  10  min 
under  conditions  typical  for  the  troposphere,  compared  to 
only  1%  when  ammonia  was  excluded.  Their  observations 
also  excluded  the  postulated  formation  of  gas-phase  prod¬ 
ucts  that  would  later  deposit  on  surfaces,  particularly  the  in¬ 
termediate  formation  of  ammonia-sulfur  dioxide  clusters  in 
the  aerosol  [89].  Catalysis  of  the  reaction  by  the  surface  of 
the  containers  used  in  the  experiments  was  similarly  elimi¬ 
nated.  However,  the  rate  of  reaction  was  dependent  on  the 
surface-to- volume  ratio  of  the  flask  and  on  the  thickness  of 
the  water  film  adsorbed  on  the  internal  surface  of  the  vessel. 

The  results  are  consistent  with  a  reaction  mechanism 
whereby  ammonia  and  sulfur  dioxide  adsorb  synergistically 
in  the  thin  film  of  water  on  the  vessel  surface  or  on  the  skin 
of  the  water  droplet.  The  high  local  concentration  of  the 
three  reactants,  NH3,  SO2,  and  O2  in  this  thin  film  results 
in  a  high  reaction  rate  and  conversion  ratio.  The  conversion 
ratio  depends  on  the  water  content  of  the  aerosol  but  not 
on  the  type  of  aerosol  (NaCl  crystal,  soot)  used  to  nucle¬ 
ate  the  cloud  drops.  Under  typical  cloud  conditions,  in  the 
presence  of  ammonia,  the  conversion  of  SO2  to  sulfate  will 
be  completed  in  less  than  5  min,  a  conclusion  confirmed  by 
actual  cloud  water  sample  measurements  [87]. 

4.2.  Reaction  of  ammonia  with  sulfur  trioxide,  SOs 

4.2.  L  Gas-phase  reaction  of  ammonia  and  SO^ 

Sulfur  trioxide  (SO3)  is  produced  in  the  atmosphere  from 
the  gas-phase  oxidation  of  SO2  by  electronically  excited 
oxygen,  by  reactive  molecules  (O3,  NO2,  N2O5),  and  by 
various  radicals  (OH,  NO3,  HO2,  RO2,  RC(0)0, . . . )  ([40] 
pp.  165-167).  Subsequent  reaction  with  atmospheric  mois¬ 
ture  produces  sulfuric  acid.  For  a  long  time  it  was  assumed 
that  the  reaction  of  SO3  with  water  should  be  so  fast  that 
reaction  with  other  atmospheric  components  could  be  ne¬ 
glected.  However,  in  1988,  Wang  et  al.  [90]  found  that  the 
gas-phase  reaction  of  SO3  and  H2O  was  much  slower  than 
generally  assumed,  with  a  second-order  rate  constant  esti¬ 
mated  to  be  5.7  X  10“^^  cm^  mol“^  s“^  The  order  of  mag¬ 
nitude  of  this  rate  constant  was  later  confirmed  by  Reiner 
and  Arnold  [91]  who  however,  reported  a  lower  value  for  the 
rate  constant:  (1.2±0.2)  x  10”^^  cm^  mol“^  s“^  Reactions 
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with  Other  atmospheric  compounds  cannot  be  dismissed  a 
priori.  In  fact,  Shen  et  al.  [92]  found  that  the  gas-phase  re¬ 
action  of  ammonia  with  SO3  was  faster  than  the  gas-phase 
reaction  of  SO3  with  water  by  almost  four  orders  of  magni¬ 
tude:  6.9x10“^^  cm^mol"^  s”^ 

As  the  temperature  decreases  (e.g.,  under  winter  condi¬ 
tions),  the  reaction  with  ammonia  becomes  competitive  with 
the  reaction  with  water.  For  example,  the  two  reactions  are 
expected  to  have  the  same  reaction  rate  at  — 10  (263 . 15  K) 
and  50%  RH. 

Under  normal  atmospheric  conditions,  the  water  concen¬ 
tration  in  the  atmosphere  is  much  higher  than  the  concen¬ 
tration  of  ammonia  and  SO3  reacts  mostly  with  water,  for 
example: 

•  At  20  °C  and  100%  humidity  the  concentration  of  water 
in  air  is  0.94  mol  m”^. 

•  The  concentration  of  ammonia  in  clean  air  is  typically  1 
ppbv  or  44.5  x  10“^  mol  m”^  and  rarely  exceeds  25  ppbv 
in  polluted  air  which  corresponds  to  1.1 1  x  10“^  mol  m”^. 

Because  the  concentration  of  water  is  at  least  6  orders  of 
magnitudes  greater  than  the  concentration  of  ammonia,  un¬ 
der  normal  atmospheric  conditions  SO3  will  react  predom¬ 
inantly  with  water,  not  with  ammonia.  However,  under  the 
conditions  expected  during  a  massive  release  of  ammonia, 
reaction  of  SO3  with  ammonia  will  be  important. 

SO3  is  a  strong  Lewis  acid  and  NH3  a  good  electron-pair 
donor.  Spectroscopic  studies  in  nitrogen  matrices  at  12  K 
have  found  evidence  for  the  formation  of  a  1:1  adduct 
through  the  donation  of  the  nitrogen  lone  electron-pair  to  the 
sulfur  on  the  SO3  moiety  [93].  More  recently,  Canagaratna 
et  al.  [94],  using  pulsed  nozzle  Fourier  transform  microwave 
spectroscopy,  were  able  to  demonstrate  the  formation  of 
this  donor-acceptor  complex  and  to  elucidate  its  structure. 

NH3(g)  +  S03(g)  ^  +H3N  -  S03-(g) 

->  H2N  -  S03H(s)  (34) 

The  N-S  bond  length  at  1.957  A  is  longer  by  0.1 86  A 
than  the  corresponding  distance  in  crystalline  sulfamic  acid 
(H2N-SO3H)  and  the  NSO  angle  at  97.6°  is  4.9°  smaller 
than  the  corresponding  angle  in  the  crystal.  These  two  ob¬ 
servations  suggest  that  the  dative  bond  is  only  partially 
formed  by  donation  of  the  lone  electron  pair  in  the  gas-phase 
adducts.  This  conclusion  is  confirmed  by  the  determination 
of  the  nitrogen  quadrupole-coupling  constant  that  indicates 
that  about  0.36  electrons  are  transferred  upon  formation  of 
the  complex.  Upon  crystallization,  the  zwitterionic  complex 
rearranges  to  give  sulfamic  acid  (Eq.  (34)).  These  results  are 
in  good  agreement  with  predictions  made  from  theoretical 
ab  initio  calculations  [95]. 

Lovejoy  and  Hanson  [96]  and  Lovejoy  [97]  have  inves¬ 
tigated  in  detail  the  kinetics  of  the  gas-phase  reaction  of 
NH3  and  SO3  in  a  nitrogen  atmosphere  over  a  wide  range  of 
temperature  (280-340  K)  and  pressure  (1.3-53  kPa),  using 
a  laminar  flow  reactor  coupled  with  a  chemical  ionization 


mass  spectrometer  for  the  detection  of  SO3  and  sulfamic 
acid.  The  second-order  rate  coefficient  for  the  NH3  and  SO3 
reaction  increases  with  increasing  concentration  of  the  car¬ 
rier  gas  N2.  The  pressure  and  temperature  dependence  of  the 
rate  coefficient  is  well  accounted  for  by  the  Troe  formalism 
developed  for  three-body  gas-phase  reactions: 


i(M)=r — ^ — ixo..? 
’  La+toOvD/fioJ 


(35) 


,3« 

where  M  refers  to  the  carrier  gas,  here  nitrogen,  ko(M)  is 
the  third  order  rate  constant  when  the  concentration  of  the 
carrier  gas  tends  toward  zero,  and  koo  is  the  second-order 
rate  constant  when  the  concentration  of  the  carrier  gas  tends 
toward  infinity.  The  experimental  data  fit  the  Troe  model  for: 

ko  (cm^  molecule”^  s“^.)  =  (3.9  ±  0.8)  x  10“^® 

koo  (cm^  molecule"^  s”^)  =  (4.7  ±  1.4)  x  10”^^ 


From  this  model,  the  effective  second-order  rate  constant  for 
the  NH3/SO3  reaction  in  1  atm  N2  at  295  K  was  calculated 
to  be  2  X  10””^^  cm^  molecule”^  s”^ 

In  the  same  study,  Lovejoy  and  co-workers  investigate 
the  fate  of  the  sulfamic  acid  formed  and  concluded  that  the 
following  mechanisms  will  control  the  loss  of  sulfamic  acid 
in  the  troposphere: 

•  Unimolecular  decomposition  yielding  back  NH3  and  SO3 : 
k^  =  9x  10’”^  s~^  for  a  pressure  of  1  atm  and  temperature 
of  281 K. 

•  Scavenging  by  aerosols. 

•  Clustering  with  sulfuric  acid  molecules  and  scavenging 
of  the  clusters  by  aerosols. 

Although  unimolecular  decomposition  dominates  under 
the  very  low  ammonia  concentrations  found  under  normal 
tropospheric  conditions,  sulfamic  acid  should  be  stable  un¬ 
der  the  high  ammonia  concentrations  found  during  massive 
release  events. 


4.2.2.  Gas-phase  reaction  of  ammonia  and  SO2  in  the 
presence  of  water 

As  mentioned  in  the  preceding  section,  when  formed  in 
the  atmosphere,  SO3  can  in  principle  react  with  either  wa¬ 
ter  or  ammonia.  The  extent  of  reaction  with  each  compound 
will  depend  on  the  relative  reaction  rates  and  concentrations 
of  the  two  species.  In  reality,  given  the  strong  physical  and 
chemical  affinities  of  these  three  compounds,  sequential  or 
concerted  series  of  reactions  among  the  three  chemicals  and 
their  intermediary  products  will  likely  determine  the  out¬ 
come  of  these  interactions. 

Tao  and  co-workers  [98,99]  used  density  functional  and 
ab  initio  molecular  orbital  calculations  to  investigate  the 
structure,  binding  energy,  and  energy  changes  associated 
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with  the  formation  and  inter-conversion  of  the  complexes 
and  products  that  can  conceivably  occur  in  the  SO3,  H2O, 
NH3  system.  Their  results  suggest  that  the  intermediate  fate 
of  SO3  in  the  atmosphere  depends  on  the  relative  concentra¬ 
tions  of  H2O  and  NH3  that  determine  whether  H2O  or  NH3 
is  initially  associated  with  SO3.  Under  normal  atmospheric 
conditions,  where  the  concentration  of  water  far  exceeds 
the  concentration  of  ammonia,  the  complex  SO3  H2O  is 
formed  initially.  In  the  presence  of  ammonia,  hydrogen 
bonded  complex  SO3  H2O  NH3  is  first  formed  that  quickly 
rearranges  to  yield  the  hydrogen-bonded  H2SO4  NH3  com¬ 
plex.  A  small  activation  energy  barrier  of  2.25kcalmol~^ 
separates  these  two  complexes.  Although  the  complex 
S03  NH3-H20  is  thermodynamically  more  stable  than  the 
complex  S03-H20-NH3,  rearrangement  does  not  occur  be¬ 
cause  of  the  large  energy  barrier  (IS.Skcalmol^^)  between 
the  two  complexes.  The  hydrogen-bonded  H2S04*NH3  com¬ 
plex  is  further  stabilized  by  association  with  two  molecules 
of  water  and  a  rapid  proton  transfer  leading  to  the  formation 
of  aqueous  ammonium  bisulfate:  HSO4”  •  NH4"'"  •  •  •  2H2O. 

In  cases  where  ammonia  is  present  in  excess  relative  to 
water,  as  would  be  the  case  for  massive  release  of  ammonia, 
the  electron  donor  acceptor  (EDA)  complex  SO3  NH3  H2O 
complex  is  formed.  This  cluster  is  thermodynamically  and 
kinetically  stable.  The  fate  of  this  complex  is  not  completely 
understood.  As  mentioned  in  the  preceding  section,  it  is 
likely  to  be  rearranged  into  sulfamic  acid  that  can  either 
dimerize  and  be  scavenged  by  aerosols,  or  form  associa¬ 
tions  with  sulfuric  acid  or  ammonium  bisulfate  aerosols,  if 
present. 

4.3,  Reaction  of  ammonia  with  sulfuric  acid  and  sulfate 
aerosols 

The  existence  of  acidic  sulfate  aerosols  in  the  atmosphere 
has  been  widely  reported  (see  Harrison  and  Kitto  [100]  for 
review  of  relevant  literature).  These  aerosols  play  an  im¬ 
portant  role  in  influencing  the  earth’s  climate  by  directly 
scattering  radiation,  and  indirectly,  by  serving  as  cloud  con¬ 
densation  nuclei  (CCN).  Because  of  their  importance,  the 
mechanism  of  formation  and  growth  of  these  aerosol  par¬ 
ticles  has  been  extensively  studied  but  is  still  incompletely 
understood. 

4.3.1.  Kinetics  of  ammonia  reaction  with  sulfuric  acid 
aerosols  under  controlled  laboratory  conditions 

From  a  phenomenological  standpoint,  the  interaction  of 
ammonia  with  sulfuric  acid-water  aerosol  is  very  similar  to 
the  interaction  of  ammonia  with  water  described  in  Section 
2.3.1  The  following  steps  are  involved: 

•  Diffusion  of  ammonia  in  the  gas  phase  to  the  liquid  sur¬ 
face. 

•  Impingement  of  the  ammonia  molecule  with  the  liquid 
surface  resulting  either  in: 


o  capture  and  adsorption  of  the  molecule  at  the  surface, 
or 

o  inelastic  scattering  away  from  the  surface. 

If  adsorption  takes  place,  the  adsorbed  molecule  can  ei¬ 
ther: 

•  desorb  firom  the  surface  and  return  to  the  gas  phase; 

•  remain  adsorbed  as  a  complex  formed  with  solvent 
molecules; 

•  react  with  chemicals  present  at  the  interface  to  form  stable 
products. 

The  products  formed  at  the  interface  can  then  migrate  by 
diffusion  into  the  bulk  of  the  liquid  opening  new  sites  for 
additional  capture  of  gas-phase  molecules. 

Consequently,  the  kinetics  of  the  overall  process  can  be 
controlled  by  one  of  the  following  steps: 

•  The  diffusion  of  the  ammonia  molecules  in  the  gas  phase. 

•  The  reaction  rate  at  the  interface. 

•  The  diffusion  of  the  products  of  reaction  away  from  the 
surface  and  into  the  bulk  of  the  liquid. 

Depending  on  the  experimental  conditions  selected,  each 
one  of  these  three  steps  can  become  the  rate-controlling  step, 
a  situation  that  was  not  initially  anticipated  and  that  caused 
a  lot  of  confusion  in  the  literature,  for  example: 

•  Robbins  and  Cadle  [101]  using  low  sulfuric  acid  con¬ 
centration  (19.7|xmolm“^),  stoichiometric  ratios  for 
H2SO4/NH3  ranging  from  0.1  to  0.48  and  drop  diam¬ 
eter  ranging  from  0.2  to  0.9  |xm,  found  that  the  initial 
reaction  was  of  first-order  with  respect  to  NH3  concen¬ 
tration  and  drop  size  but  that,  following  a  few  seconds, 
the  reaction  kinetics  did  not  follow  any  particular  or¬ 
der.  The  rate  constant  of  the  initial  reaction  was  found 
to  be  7.6  X  lO'^cm^mol"'^  s“^  (standard  deviation  of 
0.92  X  10^)  with  a  low  activation  energy  of  3  kcal  mol“^ 
Because  increasing  ammonia  partial  pressure  or  changing 
the  carrier  gas  from  nitrogen  to  helium  had  no  significant 
impact  on  the  reaction  rate,  they  concluded  that  gas-phase 
di^sion  was  not  controlling  and  that  the  overall  process 
was  controlled  by  the  rate  of  diffusion  of  the  product 
of  reaction  into  the  bulk  of  the  droplet.  They  calculated 
a  collision  efficiency  of  0.1,  i.e.,  1  out  of  10  ammonia 
molecules  hitting  the  droplet  surface  was  captured  and 
reacted  at  the  interface. 

•  The  same  authors  [102]  re-examining  earlier  kinetic  data 
developed  by  Johnstone  and  Williams  [103]  who  used 
larger  droplets  (2-5  mm),  higher  ammonia  concentrations 
(1-5%,  v/v),  and  droplets  with  acid  concentration  of  up 
to  12%  concluded  that  the  overall  reaction  rate  was  con¬ 
trolled  by  gas-phase  diffusion  and  that  every  collision  of 
ammonia  with  the  surface  resulted  in  reaction. 

All  subsequent  laboratory  investigations  [100,104-108] 
reinforce  the  importance  of  initial  conditions  in  determining 
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the  course  of  the  subsequent  reactions.  The  following  general 
conclusions  can  be  drawn: 

•  The  initial  reaction  occurs  at  the  gas-phase/liquid  inter¬ 
face  and  is  of  first  order  with  respect  to  both  gas-phase 
ammonia  concentration  and  droplet  surface  area. 

•  If  the  partial  pressure  of  NH3  is  low,  subsequent  reactions 
are  controlled  by  gas-phase  diffusion. 

•  If  the  partial  pressure  of  NH3  is  high,  subsequent  reaction 
rate  is  controlled  by  the  rate  of  diffusion  of  the  products 
of  reaction  into  the  bulk  of  the  droplet. 

•  The  uptake  coefficient  increases  as  a  function  of  acid  con¬ 
centration  and  reaches  unity  at  about  55  wt.%  H2SO4  in 
the  droplet. 

•  The  nature  of  the  product  formed  is  a  function  of  the 
stoichiometric  ratio  of  the  two  reactants. 

In  the  case  of  a  massive  release  of  ammonia  in  the  atmo¬ 
sphere,  the  very  high  concentration  of  ammonia  and  over¬ 
whelming  stoichiometric  ratio  of  ammonia  over  sulfuric  acid 
will  ensure  a  very  rapid  reaction  rate  with  sulfuric  acid 
aerosols  and  formation  of  ammonium  sulfate  (NH4)2S04  as 
the  only  product. 

43,2,  Reaction  of  ammonia  with  sulfuric  and  sulfate 
aerosols  under  normal  atmospheric  conditions 

Importance  of  this  reaction:  New  particle  formation  by 
nucleation  of  gas-phase  species  significantly  influences  the 
size  and  number  of  tropospheric  aerosols.  These  aerosols 
can  affect  the  earth’s  radiation  budget  directly  by  scattering 
solar  radiation  or  by  indirectly  serving  as  cloud  condensation 
nuclei  (CCN).  Recently,  McMurry  et  al.  [109]  and  Coffman 
and  Hegg  [110]  observed  that  the  rate  of  formation  of  new 
nanoparticles  (3-10  nm)  in  the  troposphere  of  clean,  remote 
background  regions  was  several  orders  of  magnitude  faster 
than  predicted  by  the  binary  theory  for  nucleation  of  sulfu¬ 
ric  acid  and  water.  They  suggested  that  a  ternary  nucleation 
model  involving  atmospheric  ammonia,  H2S04*NH3*H20, 
could  account  for  the  fast  initial  nucleation  rate.  Measure¬ 
ments  by  Weber  et  al.  [1 1 1]  at  a  remote  marine  site  (Mauna 
Loa  Observatory,  Hawaii)  and  at  a  remote  continental  site 
(Idaho  Hill,  Colorado),  confirmed  the  inadequacy  of  the  bi¬ 
nary  model  to  account  for  the  rate  of  nanoparticle  formation 
and  for  their  formation  at  lower  sulfuric  acid  concentra¬ 
tions  than  predicted  by  the  binary  model.  During  the  First 
Aerosol  Characterization  Experiment  (ACEl)  conducted 
in  the  sub-Antarctic  island  of  Macquarie,  south  of  Aus¬ 
tralia,  Weber  et  al.  [1 12]  measured  a  large  number  of  newly 
formed  nanoparticles  directly  downwind  from  penguin 
colonies  (source  of  ammonia).  Furthermore,  the  nanopar¬ 
ticle  concentration  was  also  correlated  with  the  gas-phase 
concentration  of  sulfuric  acid,  firmly  establishing  the  partic¬ 
ipation  of  both  sulfuric  acid  and  ammonia  in  the  formation 
of  these  ultra  small  particles.  Kerminen  et  al.  [113]  further 
showed  that  particle  growth  is  driven  initially  by  the  flux 
of  sulfuric  acid,  followed  by  thermodynamic  equilibration 
by  the  more  abundant  water  and  ammonia  vapors.  The  ini¬ 


tial  fast  formation  of  nanoparticles  is  followed  by  a  slow 
growth  phase  during  which  particles  begin  to  accumulate 
other  atmospheric  contaminants  such  as  HCl  and  HNO3 
together  with  ammonia,  resulting  in  exponential  growth 
to  particles  greater  than  50-100  nm  that  may  act  both  as 
scattering  material  and  as  cloud  condensation  nuclei. 

Apparent  reaction  rates  in  the  atmosphere:  Given  the  low 
concentration  of  ammonia  typically  present  in  the  atmo¬ 
sphere  and  given  the  potential  presence  of  several  acidic 
components  in  the  acidic  aerosols  (H2SO4,  HCl,  HNO),  a  di¬ 
rect  measurement  of  the  reaction  rate  of  ammonia  with  sulfu¬ 
ric  acid  in  the  atmosphere  is  very  difficult.  Even  so,  Harrison 
and  Kitto  [100]  carried  out  this  determination  by  measuring 
the  conversion  of  NH3  to  NH4'^  as  the  air  mass  traveled  be¬ 
tween  three  successive  distant  sampling  points  in  rural  Eng¬ 
land.  They  established  that  the  sum  of  the  concentrations 
[H++NHJ]  as  well  as  the  ratio  [H++NHJ]/[S0^"+N0^] 
remained  constant  indicating  no  significant  ammonia  input 
during  travel.  The  time  of  travel  between  locations  was  too 
short  (1-2  h)  to  permit  significant  formation  of  H2SO4  from 
SO2  oxidation  during  travel.  The  dominant  acidic  species 
in  the  aerosol  was  NH4  HSO4.  They  also  found  a  pseudo 
first-order  reaction  rate  constant  for  the  conversion  of  am¬ 
monia  to  ammonium  of  (1.5  d:  1.3)  x  10““^  s“^  between  the 
two  first  sampling  points  and  of  (1.0  ±  1.1)  x  10”"^  s”^  be¬ 
tween  the  second  and  third  sampling  points.  The  decrease 
in  the  apparent  rate  constant  as  the  reaction  progresses  is 
consistent  with  laboratory  observation  made  by  Huntzicker 
et  al.  [105]  and  McMurry  et  al.  [106]  and  was  attributed 
to  the  increased  resistance  to  neutralization  related  to  mass 
transfer  resistance  of  the  products  through  aerosol  droplets 
because  of  an  approaching  saturation  or  even  crystalliza¬ 
tion  of  the  product.  The  apparent  reaction  rate  is  positively 
correlated  with  the  acidity  of  the  aerosol,  expressed  by  the 
[H+]/[NH4+]  mol  ratio  (Eq.  (37)). 

is:  =  23x  10-5  (r  =  0.65)  (37) 

These  rate  constants  are  in  general  agreement  with  other 
determinations  of  the  ammonia  neutralization  rates  in  the 
atmosphere.  Erisman  et  al.  [114]  reported  an  average  con¬ 
version  rate  constant  of  1.0  x  lO”"^  s“^  during  daytime  pe¬ 
riods  and  5  x  10“^  s"”^  during  nighttime  periods  with  nitric 
acid  being  the  main  acidic  component  in  the  atmosphere. 
To  make  their  long-range  transport  model  fit  actual  atmo¬ 
spheric  measurements,  Asman  and  van  Jaarsveld  [115]  used 
a  value  for  ^  of  8  x  10“^  s*"^ . 

•  Products  of  reaction-aerosol  growth:  The  gas-phase  re¬ 
action  of  ammonia  and  sulfuric  acid  can  lead  to  three 
well-defined  compounds  depending  on  the  molar  ratio  of 
the  two  reactants: 

o  NH3/H2S04  =  1;  ammonium  bisulfate  NH4HSO4 
o  NH3/H2SO4  =  1.5;  letovicite  (NH4)3H(S04)2 
o  NH3/H2S04  =  2;  ammonium  sulfate  (NH4)2S04 
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These  compounds  are  stable  under  atmospheric  condi¬ 
tions.  Scott  and  Cattell  [116]  investigated  the  stability  of 
these  ammonium  sulfates  between  45  and  180  °C  using  the 
ring  oven  technique.  The  behavior  is  a  complex  function 
of  temperature  and  composition,  specifically  of  the  propor¬ 
tion  of  acidity  in  the  sample,  i.e.,  NH4HS04/(NH4)2S04 
ratio.  The  partial  pressure  of  NH3  over  the  solid  is  lowered 
in  proportion  to  its  acidity.  For  example,  at  25  °C,  a  sam¬ 
ple  containing  99.99%  (NH4)2S04  and  0.01%  NH4HSO4 
exhibits  an  NH3  partial  vapor  pressure  of  4  x  10“^^  at¬ 
mospheres  while  a  sample  containing  99%  (NH4)2S04 
and  1%  NH4HSO4  has  an  NH3  partial  vapor  pressure  of 
only  4  X  10“^^  atmospheres.  This  low  vapor  pressure  in¬ 
dicates  ammonium  sulfate,  (NH4)2S04,  will  be  the  most 
stable  form  of  solid  in  the  atmosphere,  except  for  in  the 
complete  absence  of  ammonia,  and  will  be  the  initial  sta¬ 
ble  product  formed  by  the  gas-phase  reaction  of  ammo¬ 
nia  with  sulfuric  acid  aerosols.  During  the  growth  of  the 
aerosol  particle  by  capture  of  additional  ammonia,  sulfuric 
acid,  or  water,  equilibrium  conditions  will  be  established 
that  will  change  the  initial  speciation  of  the  condensation 
nuclei. 

Depending  on  the  value  of  the  molecular  ratio,  x,  of  am¬ 
monia  and  sulfuric  acid  and  depending  on  the  trend  of  the 
humidity  in  the  atmosphere,  the  particle  of  anunonium  sul¬ 
fate  can  be  solid,  liquid,  or  it  can  be  solid-and-liquid.  The 
phase  diagram  of  the  system,  (NH4)2S04“H2S04-H20  is 
now  well  understood  as  shown  in  Fig.  2,  as  a  result  of  sev¬ 
eral  studies  aimed  at  modeling  the  formation  and  growth  of 
acidic  aerosols  [117-124]. 


Fig.  2.  Solubility  diagram  water-ammonium  sulfate-sulfuric  acid 
[117-119]. 


Depending  upon  the  beginning  composition,  mixed  sul¬ 
fate  aerosols  undergo  a  phase  transition  at  two  distinct  rel¬ 
ative  humidities  that  signal  the  onset  of  deliquescence  for 
the  two  stoichiometric  mixed  crystals,  ammonium  bisulfate, 
NH4HSO4,  and  letovicite,  (NH4)3H(S04)2.  The  growth  of 
mixed  sulfate  aerosol  particles  is  best  described  by  follow¬ 
ing  the  fate  of  a  particle  as  a  function  of  both  increasing 
X  =  NH3/H2SO4  and  increasing  relative  humidity. 

•  For  very  acidic  aerosol  particles,  jc  <  1,  and  low  humidity, 
RH  <  40%,  the  particles  consist  of  a  liquid  droplet  with 
an  increased  concentration  of  crystalline  NH4HSO4  as  x 
approaches  1 .0,  crystals  appear  when  x  reaches  about  0.2 
and  relative  humidity  exceeds  20%. 

•  For  X  =  1  a  single  stage  phase  transition  occurs  at  39.5% 
RH  and  the  crystalline  NH4HSO4  particles  grow  to  be  del¬ 
iquescent  and  commence  to  grow  by  water  condensation 
until  the  crystal  is  fully  dissolved  in  the  water  drop. 

•  For  1  <  X  <  1.5,  the  droplet  consists  of  a  solid  core  of 
a  mixed  bisulfate  and  letovicite  for  humidity’s  below  the 
deliquescent  point  for  ammonium  bisulfate,  i.e.,  39.5% 
RH.  At  relative  humidity  above  39.5%  RH,  the  parti¬ 
cle  grows  by  water  condensation,  dissolving  the  bisulfate 
crystals  and  resulting  in  an  aqueous  solution  of  bisulfate 
in  equilibrium  with  the  letovicite  crystals.  When  the  RH 
reaches  69.3%,  the  deliquescence  point  of  letovicite,  a 
second  phase  transition  begins  to  take  place  and  the  par¬ 
ticle  grows  by  water  condensation  until  all  the  letovicite 
is  dissolved. 

•  For  1.5  <  x  <  2,  the  particle  consists  of  a  crystalline  core 
of  mixed  letovicite  and  sulfate  until  the  RH  reaches  the 
deliquescent  point  of  letovicite,  69.3%,  where  the  particle 
grows  by  water  condensation  and  the  letovicite  dissolves 
leaving  ammonium  sulfate  crystals  in  equilibrium  with 
the  solution.  Ammonium  sulfate  deliquesces  at  about  80% 
RH. 

In  the  case  of  a  massive  release  of  ammonia  in  the  at¬ 
mosphere,  the  stoichiometric  ratio,  x,  will  be  significantly 
greater  than  two.  The  initial  product  of  the  gas-phase  reac¬ 
tion  of  H2SO4  aerosols  and  NH3  will  be  ammonium  sulfate 
together  with  sulfamic  acid,  at  a  very  low  humidity  level. 
For  atmospheric  relative  humidity  below  80%,  nanoparti¬ 
cles  of  ammonium  sulfate  should  be  formed  from  reaction 
with  pre-existing  sulfuric  acid  aerosols  together  with  some 
sulfamic  acid  produced  by  the  reaction  of  ammonia  with 
freshly  formed  SO3  resulting  from  the  ammonia-catalyzed 
oxidation  of  SO2.  For  RH  above  80%,  the  particles  will  grow 
by  water  condensation  until  thermodynamic  equilibrium  is 
reached  for  the  system  [(NH4)2S04]  solid-saturated  aqueous 
solution-NH3  gas.  Because  of  the  extraordinary  conditions 
that  it  represents,  this  system  has  not  yet  been  investigated. 
It  is  conceivable,  however,  that  these  particles  could  grow 
to  the  50-100  nm  size  characteristics  of  cloud  condensation 
nuclei  by  water  condensation  and/or  by  accretion  of  other 
products  of  ammonia  reaction  with  atmospheric  pollutants, 
for  instance,  ammonium  nitrate  and/or  ammonium  chloride. 


46 


JJ.  Renard  et  al/ Journal  of  Hazardous  Materials  BIOS  (2004)  29-60 


5.  Reaction  of  ammonia  with  atmospheric  nitrogen 
oxides  and  nitric  acid 

Nitric  oxide,  NO,  and  nitrogen  dioxide,  NO2,  are  primary 
pollutants  emitted  by  both  mobile  combustion  sources  and 
emitted  by  stationary  combustion  sources.  Under  normal  at¬ 
mospheric  conditions,  NO  is  converted  into  NO2  within  a 
small  number  of  hours.  Through  a  variety  of  mechanisms  in¬ 
volving  both  oxidation  and  hydrolysis,  NO2  rapidly  ends  up 
as  either  nitric  acid  or  as  peroxyacetyl  nitrate  (PAN)  and  to 
a  slighter  extent  may  end  up  as  a  variety  of  organic  nitrates 
([39]  Chapter  8  and  Section  11-B). Ammonia  is  known  to 
reduce  nitrogen  oxides  back  to  nitrogen.  However,  this  re¬ 
action  requires  high  temperatures  and/or  catalysts  to  occur 
at  measurable  rates  and  its  occurrence  under  atmospheric 
conditions  is  unlikely  and  in  fact  has  never  been  reported. 
To  our  knowledge,  the  reaction  of  ammonia  with  PAN  has 
not  yet  been  studied.  It  is  conceivable,  however,  that  these 
two  chemicals  would  react,  particularly  at  the  surface  of  wa¬ 
ter  droplets  according  to  the  well-known  alkaline  hydrolysis 
reaction  of  PAN  ([39],  pp.  541-542)  (Eq.  (38)). 

CH3C(=0)0N02  -f  2NH3  +  H2O 

^  NH4CH3CO2  +  NH4NO2  +  5O2  (38) 

To  our  knowledge,  the  presence  of  ammonium  nitrite, 
NH4NO2,  or  ammonium  acetate,  NH4CH3CO2,  has  not 
been  reported  in  atmospheric  aerosols.  However,  normal 
atmospheric  conditions  do  not  support  simultaneous  ele¬ 
vated  concentrations  of  both  PAN  and  of  ammonia  in  the 
same  environment.  PAN  tends  to  be  predominating  in  ur¬ 
ban  environments,  and  ammonia  tends  to  be  predominating 
in  rural  environments.  In  the  case  of  a  massive  release  of 
ammonia  in  an  urban  environment  a  reaction  of  PAN  with 
ammonia  may  possibly  occur.  The  average  concentrations 
of  PAN  range  from  2  ppbv  in  rural  environments  to  70  ppbv 
in  heavily  polluted  urban  environments  ([39],  p.  369). 

Aside  from  PAN,  nitric  acid  is  the  final  product  for  all 
nitrogen  oxides  formed  in  the  atmosphere.  Nitric  acid  has 
a  great  deal  higher  vapor  pressure  than  does  sulfuric  acid 
under  typical  atmospheric  conditions.  It  is  therefore  present 
in  significant  concentrations  in  the  gas-phase  rather  than 
being  present  in  an  aerosol  form,  as  is  the  case  for  sulfuric 
acid. 

5.7.  Gas-phase  reaction  of  ammonia  and  nitric  acid 

The  reaction  of  HNO3  and  NH3  is  a  prototypical  acid-base 
neutralization,  i.e.,  proton  transfer  reaction,  which,  for  all 
practical  purposes,  occurs  instantaneously  in  aqueous  solu¬ 
tions.  However,  there  is  no  experimental  evidence  that  this 
reaction  takes  place  in  the  gas  phase.  Tao  and  co-workers 
[125,126]  have  used  high  level  ab  initio  calculations  to  inves¬ 
tigate  the  equilibrium  structures,  binding  energies,  and  har¬ 
monic  frequencies  of  the  potential  ternary  clusters  formed 
by  interaction  of  nitric  acid,  ammonia,  and  water  molecules 


in  the  gas  phase.  They  mapped  the  potential  energy  surfaces 
along  the  proton  transfer  pathway  to  identify  the  preferred 
reaction  path. 

Their  results  indicate  that  in  the  absence  of  water,  the  ni¬ 
tric  acid-ammonia  system  exists  as  hydrogen  bonded,  with 
nitric  acid  acting  as  the  hydrogen  bond  donor  and  ammonia 
as  the  acceptor.  Addition  of  one  water  molecule  stabilizes 
the  hydrogen-bonded  complex  but  is  not  enough  to  stabilize 
the  ion-pair  that  would  result  from  proton  transfer.  The  sit¬ 
uations  change  with  the  addition  of  a  second  molecule  of 
water  that  produces  additional  stabilization  energy,  which  in 
turn  helps  stabilize  the  ion-pair  and  thus  induces  the  trans¬ 
fer  of  a  proton  from  the  nitric  acid  to  the  ammonia  to  form 
the  ammonium  ion.  A  third  water  molecule  contributes  to 
further  stabilization  of  the  ion-pair. 

Under  normal  atmospheric  conditions,  the  concentration 
of  water  in  the  air  is  several  orders  of  magnitudes  higher  than 
concentrations  of  either  nitric  acid  or  ammonia,  therefore; 
one  can  expect  the  formation  of  water-stabilized  ammonium 
nitrate  clusters  in  the  atmosphere. 

5.2.  Atmospheric  reaction  of  ammonia  and  nitric  acid 

The  reaction  of  ammonia  with  nitric  acid  yields  ammo¬ 
nium  nitrate  that  is  in  equilibrium  with  its  gaseous  con¬ 
stituents,  whether  ammonium  nitrate  is  a  crystalline  solid  or 
a  solute  in  an  aqueous  solution  (Eq.  (39)). 

NH4N03(s or  aq)  ^  NH3(g)  +  HN03(g)  (39) 

The  stability  of  ammonium  nitrate  has  been  the  object 
of  several  studies  [127,128].  NH4NO3  does  not  start  its 
well-known  rapid,  self-accelerating  decomposition  to  N2O 
and  water  until  heated  to  260-280  °C.  Below  this  tempera¬ 
ture  it  undergoes  reversible  decomposition  into  its  compo¬ 
nents  as  shown  in  Eq.  (39),  whether  in  the  form  of  melted 
salt  above  its  melting  point  of  169.6  °C  or  in  the  form  of  a 
crystalline  solid  or  aqueous  solution  at  lower  temperatures. 
Brandner  et  al.  [129]  extended  the  temperature  range  over 
which  the  vapor  pressure  above  ammonium  nitrate  was  mea¬ 
sured  and  calculated  heat  of  vaporization  of  42.7  kcal  mol”^ 
for  the  solid  and  39.9  kcal  mol*"  ^  for  the  liquid,  in  reason¬ 
able  agreement  with  the  theoretical  value  of  42.1  kcal  mol”^ 
computed  from  thermodynamic  data. 

By  the  mid  1970s,  a  number  of  sampling  studies  had 
shown  that  atmospheric  aerosols  often  contained  both  am¬ 
monia  and  nitric  acid.  Stelson  et  al.  [130]  reviewed  and 
analyzed  these  results  and  demonstrated  that  ammonia  and 
nitric  acid  were  found  in  a  roughly  1:1  mole  ratio  and  the 
product  of  their  gas-phase  concentration  was,  within  experi¬ 
mental  error,  equal  to  the  expected  equilibrium  constant  for 
Eq.  (39).  This  observation  proved  that  ammonium  nitrate  ex¬ 
isted  as  aerosol  particulates  in  the  atmosphere  where  it  was 
in  equilibrium  with  the  vapor  pressure  of  its  constituents  and 
it  was  not  formed  as  an  artifact  on  the  filter  media  used  to 
collect  the  particulates.  These  conclusions  were  confirmed 
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Table  13 


Equilibrium  constant  for  solid  ammonium  nitrate 

Reference 

Equation  (Kc  in  ppm^) 

Kcat  273.15  K 

Kcat  298. 13  K 

Kcat  303.15  K 

Brandner  et  al.  [129] 

Stelson  et  al.  [130] 

Stelson  and  Seinfeld  [135] 

InKc  =62.296  -  (21.510/7) 

In  Kc  =70.62  -  (24.090/7)  -  6.041n(77298.15) 
InKc  =70.78  -  (24.220/7)  -  6.101n(r/298.15) 

7.16  X  10-® 

1.37  X  10-® 

9.97  X  10“* 

5.28  X  10"’ 

3.80  X  10-5 

2.89  X  10-5 

1.73  X  10-^ 
1.30  X  10-'' 

9.9  X  10-5 

by  additional  measurements  by  Doyle  et  al.  [131]  using 
kilometer-length  Fourier-transform  infra  red  spectroscopy. 

Nevertheless,  significant  problems  continued  to  be  re¬ 
ported  with  the  accuracy  and  reliability  of  atmospheric 
aerosol  nitrate  concentrations.  For  example,  Appel  et  al. 
[132]  found  the  atmospheric  concentrations  of  ammonia 
and  nitric  acid  at  Pittsburgh  (CA)  yielded  values  generally 
below  those  needed  for  agreement  with  the  equilibrium 
constant  of  Eq.  (39).  Forrest  et  al.  [133]  demonstrated  the 
NH4NO3  deposited  on  filter  media  was  lost  to  volatilization 
for  relative  moisture  below  60%  but  no  measurable  loss 
occurred  at  100%  relative  humidity.  To  gain  an  improved 
understanding  of  the  atmospheric  behavior  of  anunonium 
nitrate  Stelson  and  Seinfeld  undertook  extensive  studies 
of  the  factors  determining  the  equilibrium  constant  of 
NH4NO3  solid  [134]  and  in  acidic  aqueous  solutions  [135]. 

5.2.7.  Influence  of  temperature  on  solid  ammonium  nitrate 
decomposition 

As  shown  in  Table  13,  the  equilibrium  constant,  Kc,  for 
the  decomposition  of  solid  ammonium  nitrate  into  its  con¬ 
stituents  is  a  strong  function  of  temperature.  The  equation 
developed  by  Stelson  and  Seinfeld  [135]  uses  the  most  recent 
values  for  thermodynamic  data  and  takes  the  heat  capaci¬ 
ties  of  the  chemical  species  into  account.  It  should  therefore 
give  the  most  accurate  values  for  Kc. 

5.2.2.  Influence  of  relative  humidity  on  the  decomposition 
of  ammonium  nitrate 

Stelson  and  Seinfeld  [134]  have  critically  reviewed  the 
published  experimental  data  on  the  relative  humidity  of  del¬ 
iquescence  for  ammonium  nitrate.  In  general,  most  authors 
agree  that  at  25  °C,  NH4NO3  deliquesces  at  a  relative  hu¬ 
midity  of  62%.  Reports  for  start  of  deliquescence  at  relative 
humidity  as  low  as  30%  [136]  can  probably  be  attributed 
to  the  presence  of  impurities,  particularly  to  the  presence  of 
sodium  nitrate. 

The  point  of  deliquescence  is  a  strong  function  of  tem¬ 
perature.  Dingermans  [137]  found  that  experimental  data  for 
the  deliquescence  of  NH4NO3  fit  the  least  square  expression 
shown  in  Eq.  (40): 

856.23  ±  13.25 

ln(RHD)  = - - - -i-  1.2306  ±  0.0439, 

RHD298.15  K  =  60.49%  (40) 

The  thermodynamic  derivation  of  the  temperature  depen¬ 
dence  of  the  relative  humidity  at  deliquescence  yields  a  very 
similar  relationship  (Eq.  (41)). 


ln(RHD)  =  +  1.7037,  RHD298.15  k  =  62.24% 

(41) 

Eq.  (41)  being  on  a  strong  theoretical  footing  and  giving  a 
result  closer  to  the  best  defined  experimental  value  in  the 
normal  range  of  atmospheric  temperatures  should  probably 
be  preferred. 


5.2.5.  Aqueous  solubility  limit  of  ammonium  nitrate 
The  solubility  limit  of  NH4NO3  in  aqueous  solution  is  also 
an  inverse  function  of  temperature.  The  results  of  Stephen 
and  Stephen  [138]  give  the  least  square  regression,  Eq.  (42), 
for  the  dependence  of  the  molality  of  the  saturated  aqueous 
solution  as  a  function  of  temperature: 


ln(/n)  = 


-1837.3  ±18.00 


+  9.4235  ±  0.0602, 


"1298.15  K  =  26.08  mol  kg  ^  water 


(42) 


The  calculated  value  at  298.15  K  agrees  very  well  with  the 
corresponding  experimental  value  of  Hamer  and  Wu  [139]: 
25.954 mol kg~^  water. 

Stelson  and  Seinfeld’s  thermodynamic  derivation,  using 
the  Hamer  and  Wu  value  for  determining  the  integration 
constant,  yields  Eq.  (43)  that  is  very  close  to  the  experimental 
Eq.  (42). 

-1600 

ln(m)  =  — ^  ±  8.6228  (43) 


5.2.4.  Ammonium-nitric  acid  equilibrium  product  over 
aqueous  ammonium  nitrate  solutions 

Stelson  and  Seinfeld  [135]  have  made  an  extensive  ther¬ 
modynamic  analysis  of  the  NH3-NO3H-NH4NO3-H2O 
system.  Their  conclusions  can  be  summarized  as  follows: 

•  It  is  essential  to  consider  the  non-ideality  of  the  concen¬ 
trated  solutions  near  the  point  of  deliquescence.  Failure  to 
do  so  leads  to  an  order  of  magnitude  error  on  the  calcu¬ 
lation  of  the  equilibrium  product  and  of  as  much  as  20% 
error  on  the  calculation  of  the  point  of  deliquescence. 

•  The  product  of  ammonia  and  nitric  partial  pressure  above 
the  solution  is  inversely  related  to  relative  humidity  but  is 
virtually  independent  of  pH  for  solutions  with  pH  between 
1  and  7.  Below  pH  1,  the  influence  of  undissociated  nitric 
acid  must  be  taken  into  account.  At  pH  higher  than  7, 
dissolved  ammonia  complicates  the  theoretical  treatment 
of  the  solution. 

•  Unreactive  solutes  such  as  ammonium  sulfate  (see  next 
section)  lower  the  water  vapor  pressure  but  have  no  effect 
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Fig.  3.  Equilibrium  of  ammonia  and  nitric  acid  over  ammonium  nitrate 
solutions  [135]. 


on  the  ammonia-nitric  acid  vapor  pressure  product.  It  fol¬ 
lows  that  the  ammonia-nitric  acid  vapor  pressure  product 
will  be  different  at  a  given  relative  humidity  from  what 
it  would  have  been  at  the  same  relative  humidity  in  the 
absence  of  the  additional  solute. 

The  predictive  equations  developed  by  Stelson  and  Sein¬ 
feld  are  complex  and  have  no  direct  application  for  our  pur¬ 
pose.  The  overall  behavior  of  the  NH3-NO3H-NH4NO3- 
H2O  system  is  best  illustrated  in  Fig.  3  below: 

For  a  given  temperature,  the  NH4NO3  dissociation  con¬ 
stant  is  invariant  below  the  relative  humidity  of  deliques¬ 
cence.  The  on-set  of  deliquescence  occurs  for  increasing 
relative  humidity  as  the  temperature  decreases.  For  a  given 
temperature,  above  the  relative  humidity  of  deliquescence, 
the  equilibrium  constant  for  NH4NO3  dissociation  decreases 
rapidly  particularly  as  relative  humidity  approaches  100%. 
The  discontinuity  shown  around  the  deliquescence  point  has 
several  possible  causes:  use  of  different  data  bases  for  ther¬ 
modynamic  data  for  solids  and  solutions,  errors  introduced 
by  interpolation  of  data,  derivation  of  polynomial  expres¬ 
sions  for  heat  capacities,  etc.  Nevertheless,  the  behavior 
shown  in  Fig.  3  constitutes  an  accurate  representation  of  the 
real  behavior  of  the  system.  For  example,  it  accounts  for 
the  observation  that  at  98%  humidity  and  25  °C,  the  mass 
concentration  of  NH3  plus  NO3H  in  the  gas  in  phase  equi¬ 
librium  with  an  aqueous  solution  of  ammonium  nitrate  is 
only  1.9jxgm“^  while  17.9  |xgm”^  would  be  expected  if 
ammonium  nitrate  were  present  as  a  solid. 


5.2.5.  Kinetics  of  evaporation  of  ammonia  nitrate  aerosols 

Larson  and  Taylor  [140]  were  the  first  to  investigate  the 
kinetics  of  evaporation  of  NH4NO3  from  aerosols.  In  their 
laboratory  diffusion  stripper,  they  used  very  small  droplets 
(0.4  |uum)  of  aqueous  solutions  of  NH4NO3  and  continu¬ 
ous  trapping  of  emitted  NH3  and  HNO3.  The  measured 
rates  of  size  reduction  for  the  aerosol  particles  agreed  well 
with  those  predicted  from  gas  diffusion  and  concentrated 
solution  chemistry.  They  concluded  that  the  “evaporation 
coefficient”  for  both  NH3  and  HNO3  was  unity,  implying 
negligible  resistance  for  the  transport  of  molecules  across 
the  vapor-liquid  interface. 

A  very  different  picture  emerges  from  results  obtained 
by  later  investigators  particularly  when  actual  atmospheric 
data  were  collected.  Richardson  and  Hightower  [141]  levi¬ 
tated  larger  particles  (2.5-3  [jum)  of  either  solid  or  aqueous 
solutions  of  ammonium  nitrate  in  a  vacuum  in  an  electric 
quadrupole  trap  and  observed  the  rate  of  changes  of  the 
particles’  mass  and  size.  They  observed  that  fresh  particles 
evaporate  at  a  rate  of  —0.23  x  10“^  |xms“^  but  that  this  rate 
decreases  to  —0.06  x  10^"^  fims“^  after  4  h.  For  these  data 
to  agree  with  thermodynamic  predictions,  mass  accommo¬ 
dation  coefficients  (evaporation  coefficients)  of  0.02  must  be 
used  for  the  initial  evaporation  period,  decreasing  to  0.004 
after  4  h.  This  would  imply  that  that  transport  between  solid 
and  vapor  phases  is  significantly  inhibited. 

Attempts  to  build  numerical  models  to  account  for  the 
measured  concentrations  of  NH3,  HNO3,  and  NH4NO3  in 
the  atmosphere  have  been  largely  unsuccessful  [142]  as  were 
attempts  to  identify  the  rate-determining  step  in  the  mecha¬ 
nisms  involved  [143].  Nonetheless,  a  general  understanding 
of  the  factors  involved  emerges  from  several  other  studies 
[144-147]  allowing  a  qualitative  explanation  of  the  trends 
observed.  Deviations  from  predicted  values  based  on  ther¬ 
modynamics  occur  when  not  enough  time  is  available  for 
the  system  to  reach  equilibrium  after  a  change  in  the  con¬ 
centration  of  one  of  the  chemicals  involved,  or  after  a  rapid 
change  in  atmospheric  conditions  particularly  temperature 
and/or  relative  humidity.  Changes  of  concentrations  can  oc¬ 
cur  because  of  the  proximity  of  a  source  of  ammonia  (cattle 
feedlot,  chicken  or  pig  farm)  whose  contribution  may  vary 
depending  on  wind  direction  [144].  Another  example  is  the 
observed  vertical  gradient  of  nitric  acid  and  nitrate  aerosols. 
Near  the  ground,  the  rapid  dry  deposition  of  nitric  acid  de¬ 
pletes  the  concentration  of  HNO3,  triggering  the  decompo¬ 
sition  of  ammonium  nitrate  particles,  an  effect  not  observed 
at  higher  altitudes.  The  situation  is  further  complicated  by 
the  presence  of  other  anions,  particularly  sulfate  and  chlo¬ 
rides  that  compete  with  nitrate  for  ammonia.  Although  this 
general  framework  of  interactions  provides  a  qualitative  ex¬ 
planation  for  actual  atmospheric  observations,  the  nature  of 
the  kinetic  constraints  preventing  rapid  equilibration  of  the 
system  remains  unknown.  However,  in  the  case  of  a  mas¬ 
sive  release  of  ammonia  in  the  atmosphere,  the  local  high 
concentration  of  ammonia  is  expected  to  overcome  any  ki¬ 
netic  limitations  and  to  drive  the  system  rapidly  to  the  com- 
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plete  conversion  of  nitric  acid  to  ammonium  nitrate  whose 
decomposition  will  be  inhibited  by  the  large  excess  of  am¬ 
monia  present. 

5.i.  Interactions  of  sulfuric  acid,  nitric  acid,  and  ammonia 
in  atmospheric  aerosols 

Ammonium  sulfate  and  nitrate  are  ubiquitous  components 
of  atmospheric  aerosols,  having  been  observed  in  urban  and 
rural  environments  in  the  USA  and  abroad  (see  [148]  for 
review  of  relevant  literature).  Consequently,  a  substantial 
amount  of  effort  has  been  devoted  to  the  investigation  and 
modeling  of  the  NH4N03-(NH4)2S04“NH3-HN03-H20 
system. 

The  large  difference  in  the  stability  of  the  two  salts  is  one 
of  the  distinctive  characteristics  of  this  system: 


NH4N03(s)  ^  NH3(g)  +  HN03(g), 

K  =  3.03  X  10-*’  (44) 

NH4N03(aq)  NHsCg)  +  HN03(g), 

K  =  2.71  X  10"*’  (45) 

(NH4)2S04(s)  2NH3(g)  +H2S04(g), 

K  =  2.33  X  10"^®  (46) 

(NH4)2S04(aq)  2NH3(g)  +  H2S04(g), 

K  =  2.62  X  10"^®  (47) 


Because  of  the  very  low  dissociation  constant,  the  disso¬ 
ciation  of  ammonium  sulfate  is  negligible  under  most  atmo¬ 
spheric  conditions. 

5,3.  L  Thermodynamic  analysis  of  the 
NH^N0^/(NH^)2S0^H2/HN0^/H20  system 

Seinfeld  and  co-workers  [148-151]  conducted  a  com¬ 
prehensive  investigation  of  the  thermodynamics  of  the 
NH4N03/(NH4)2S04/NH3/HN03/H20  system,  following 
the  pioneering  works  of  Tang  [136],  Seigneur  et  al.  [152], 
and  Tanner  [153].  As  could  be  expected  from  such  a 
multi-component  system,  the  situation  is  complex,  partic¬ 
ularly  when  attempting  to  build  a  predictive  quantitative 
model.  The  two  main  determining  variables  are  the  molecu¬ 
lar  ratio  of  the  ionic  species  (NH4'^,  NO3",  and  S04^”) 
and  the  relative  humidity.  The  ratio  of  the  ionic  species  de¬ 
termines  the  chemical  entities  present,  the  relative  humidity 
determines  what  phase  they  will  be  present  in,  i.e.,  the 
phase  transition  from  a  solid  particle  to  a  liquid  drop  occurs 
when  the  relative  humidity  reaches  a  certain  critical  value 
corresponding  to  the  water  activity  of  the  saturated  solution. 
These  saturated  solutions  have  very  high  molality  values  and 
therefore  behave  in  a  strongly  non-ideal  behavior.  Chan  et  al. 
[15 1]  have  demonstrated  that  the  three  models  of  mixed  elec¬ 
trolyte  solutions,  the  Zdanovskii-Stokes-Robinson  (ZSR), 
the  Kusik  and  Meissner  (KM),  and  the  Pitzer  models,  give 


a  reasonably  accurate  description  of  the  system,  the  ZSR 
model  giving  the  most  consistent  prediction  over  the  widest 
range  of  composition  and  relative  humidity.  The  deviations 
between  the  prediction  of  the  model  and  experimental  data 
are  attributed  to  the  binary  and  ternary  solute-solute  inter¬ 
actions  that  are  not  taken  into  account  by  these  models.  In 
general,  the  maximum  deviation  does  not  exceed  6%.  Sein¬ 
feld  et  al.  used  the  KM  empirical  mixing  rule  in  their  ther¬ 
modynamic  analysis  of  the  system.  A  list  of  temperature  de¬ 
pendence  of  the  various  equilibrium  constants  and  chemical 
potentials  can  be  found  in  Basset  and  Seinfeld  [148].  The 
main  conclusion  reached  by  these  analyses,  particularly  as 
it  applies  to  high  concentrations  of  gas-phase  ammonia  rel¬ 
ative  to  sulfuric  and  nitric  acid  concentrations  (massive  re¬ 
lease  of  ammonia  in  the  atmosphere),  can  be  summarized  as 
follows: 

•  Only  NH4NO3  and  (NH4)2S04  will  be  present  in  the  liq¬ 
uid  phase. 

•  As  the  nitrate  to  sulfate  ratio  increases  from  0  to  1,  the  del¬ 
iquescence  point  decreases  from  80  to  62%  RH  at  25  °C. 
The  decrease  in  the  deliquescence  point  is  small  until  the 
molecular  ratio  of  nitrate/sulfate  reaches  1/4,  but  acceler¬ 
ates  thereafter,  reaching  about  75%  RH  for  nitrate/sulfate 
ratio  of  1/1, 70%  RH  for  3/1  until  reaching  the  final  value 
of  62%  RH  for  pure  ammonium  nitrate  as  illustrated  in 
Fig.  4.  Therefore,  for  relative  humidity  below  62%,  the 
product  of  reaction  will  be  solid  particles  of  ammonium 
nitrate,  sulfate,  or  their  double  salts,  depending  on  the 
molecular  ratio  of  nitric  and  sulfuric  acid  in  the  atmo¬ 
sphere.  Above  80%,  the  salts  will  be  in  the  form  of  a 
liquid  droplet  with  the  salt  concentration  depending  on 
the  prevailing  relative  humidity.  For  relative  humidity  be¬ 
tween  80%  and  62%  the  most  common  occurrence  will 
be  a  crystalline  core  particle  surrounded  by  a  saturated 


RELATIVE  HUMIDITY.  % 

Fig.  4.  The  Effect  of  (NH4)2S04  on  the  relative  humidity  dependence  of 
the  NH4NO3  dissociation  constant  [149]. 
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solution  of  the  salts,  the  nature  of  the  solid  and  dissolved 
salts  depending  on  the  nitrate-to-sulfate  ratio. 

•  The  phase  diagram  of  the  system  [148]  shows  the  for¬ 
mation  of  two  double  salts  2NH4NO3'-0^1H4)2SO4  and 
3NH4N03-(NH4)2S04. 

•  As  long  as  the  high  partial  pressure  of  gas-phase  ammo¬ 
nia  persists,  ammonium  nitrate  will  not  volatilize  and  the 
partial  pressure  of  nitric  acid  will  be  negligible.  When 
the  cloud  of  ammonia  dissipates,  the  ammonium  nitrate 
will  decompose  into  NH3  and  HNO3  until  the  product  of 
the  concentration  of  the  two  gases  reaches  its  equilibrium 
value  for  the  prevailing  temperature  and  relative  humidity. 

5.5.2.  Kinetics  of  particle  growth/decay  in  the  interactions 
ofNH4N03/(NH4)2S04/NH3/HN03/H20  system 
Hightower  and  Richardson  [154]  have  investigated  the 
evaporation  rate  of  single  solid  particles  with  initial  molec¬ 
ular  composition  of  80%  NH4NO3/20%  (NH4)2  SO4  using 
levitation  in  an  electric  quadrupole  trap.  When  the  evap¬ 
oration  was  conducted  in  vacuum,  the  evaporation  rate  of 
NH4NO3  is  as  much  as  360  times  higher  than  the  evapora¬ 
tion  rate  of  pure,  crystalline  NH4NO3.  This  suggests  that  the 
mixed  solid  acts  as  a  non-ideal  solution,  not  a  mixed  crys¬ 
tal,  and  that  the  effect  of  sulfate  is  to  prevent  crystallization 
thereby  causing  ion  mobility  and  mass  accommodation  co¬ 
efficient  (sticking  or  evaporation  coefficient)  to  remain  high. 
In  the  presence  of  water  vapor  below  the  deliquescence  point 
of  ammonium  nitrate,  62.2%  RH,  the  rate  of  NH4NO3  evap¬ 
oration  passes  through  a  maximum  for  a  relative  humidity 
of  about  40%  at  which  point  the  evaporation  rate  is  about 
10  times  faster  than  the  evaporation  rate  at  the  deliques¬ 
cent  point.  Although  the  particles  do  not  gain  any  measur¬ 
able  weight  when  water  vapor  is  introduced,  Hightower  and 
Richardson  propose  that  the  solid  particle  becomes  coated 
with  a  superficial  layer  of  water,  a  few  molecules  thick.  The 
ions  leave  the  solid  and  enter  the  water  layer  as  fully  sol¬ 
vated  species  whose  evaporation  rate  is  governed  by  the  ther¬ 
modynamics  of  supersaturated  solutions  rather  than  that  of 
mixed  solids.  For  40%  RH,  the  “water  coating  layer”  would 
be  optimal  for  maximum  evaporation  rate. 

In  their  study  of  the  thermodynamics  of  the  NH4NO3/ 
(NH4)2S04/NH3/HN03/H20  system,  Bassett  and  Seinfeld 
[150]  have  investigated  the  impact  of  particle  size  on  the 
equilibrium  state,  composition,  and  physical  state  of  the 
aerosol  particles.  In  particular,  they  have  considered  the  po¬ 
tential  effect  of  the  Kelvin  effect  that  governs  the  vapor 
pressure  of  species  in  a  drop  as  a  function  of  the  radius  of 
the  drop:  the  sharper  the  curvature  of  the  drop,  the  higher 
the  vapor  pressure  of  the  species  above  the  drop.  They  con¬ 
cluded  that  the  Kelvin  effect  had  a  significant  impact  on 
the  distribution  of  the  species  by  particle  size.  Ammonium 
sulfate  being  essentially  non-volatile  will  have  a  size  dis¬ 
tribution  controlled  by  gas-phase  diffusion  and  will  tend  to 
accumulate  in  small  size  particles.  Nitrates  being  substan¬ 
tially  more  volatile  than  sulfates  will  tend  to  evaporate  from 
small  particles  and  deposit  on  large  particles  where  surface 


curvature  effects  on  vapor  pressure  are  minimal.  The  net 
effect  is  a  mass  transfer  of  anunonium  nitrate  from  small 
to  large  particles  and  the  creation  of  a  bimodal  distribution 
with  very  fine  particles  where  ammonium  sulfate  predomi¬ 
nates  on  one  end  and  larger  particles  with  a  core  of  ammo¬ 
nium  sulfate  surrounded  by  deposited  ammonium  nitrate  at 
the  other  end.  This  type  of  distribution  has  been  observed 
repeatedly  when  sampling  atmospheric  aerosols  (see  refer¬ 
ence  [39]  pp.  786-790  for  review  of  relevant  data). 

As  mentioned  in  Section  5.2,  in  actual  atmospheric 
measurements,  deviation  from  equilibrium  concentrations, 
product  phase  and  composition  are  often  observed  (e.g., 
[153,155]).  These  deviations  are  mostly  due  to  the  inability 
of  the  system  to  reach  equilibrium  after  rapid  changes  of 
temperature,  relative  humidity,  or  gas-phase  composition, 
mainly  because  of  slow  diffusion  of  species  within  the  bulk 
of  the  particle  [148]. 

6.  Reaction  with  hydrochloric  acid 

The  level  of  HCl  in  the  troposphere  is  highly  variable 
and  not  well  documented  ([39],  p.  678).  The  main  source 
stems  from  the  reaction  of  sea-salt  aerosols  with  strong  acid 
aerosols  such  as  sulfuric  and  nitric  acids.  Anthropogenic 
sources  such  as  waste  incineration  and  automobile  emissions 
provide  the  balance  for  the  observed  HCl  concentration  in 
the  atmosphere.  Therefore,  HCl  should  be  expected  to  be 
prevalent  in  coastal  areas  but  the  observed  levels  seldom 
exceed  the  low  ppb  range. 

6.7.  Gas-phase  reaction  of  ammonia  and  hydrochloric  acid 

Like  nitric  and  sulfuric  acid,  hydrochloric  acid  is  a  strong 
acid  that  reacts  extremely  rapidly  with  ammonia  in  aqueous 
solutions,  yielding  ammonium  chloride. 

Early  investigators  using  measurements  of  vapor  densi¬ 
ties  and  dissociation  pressure,  as  well  as  thermodynamic 
calculations,  concluded  that  in  the  gas  phase,  ammonium 
chloride  was  completely  dissociated  into  HCl  and  NH3 
[156].  This  commonly  held  view  was  challenged  by  ab 
initio  calculation  of  the  electronic  structure  of  the  NH4CI 
molecule  predicting  a  minimum  in  the  NH4CI  potential 
energy  surface  that  would  be  19kcalmol'”^  lower  than  the 
energy  of  its  dissociation  products  [157].  Shortly  thereafter, 
Golfinger  and  Verhagen  [158]  using  quadrupole  mass  spec¬ 
troscopy  provided  unambiguous  evidence  for  the  presence 
of  undissociated  deuterated  ammonium  chloride  in  the  va¬ 
por  effusing  from  a  Knudsen  cell.  Shibata  [159]  further 
demonstrated  the  presence  of  undissociated  ammonium 
chloride  in  the  gas  phase,  using  electron  diffraction  spec¬ 
troscopy.  De  Kruif  [160]  using  a  torsion  effusion  technique 
where  electromagnetic  forces  are  used  to  compensate  for 
the  recoil  of  vapor  molecules  effusing  from  a  cell  in  which 
ammonium  chloride  is  vaporized,  determined  that  approxi¬ 
mately  85%  of  NH4CI  is  dissociating  into  NH3  and  HCl  in 
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the  gas  phase  at  352  K  (79  °C).  The  degree  of  dissociation 
increases  only  slowly  with  increasing  temperature. 

Advanced  ab  initio  calculations  by  Tao  and  co-workers 
[161-164]  support  the  formation  and  stability  of  NH4CI 
in  the  gas  phase.  The  single  NH3-HCI  complex  is  hydro¬ 
gen  bonded  only.  The  proton  transfer  does  not  come  about 
until  two  molecules  of  water  are  present  or  until  at  least 
a  cluster  of  two  NH3-HCI  complexes  form.  The  clusters 
of  two  or  four  NH3-HCI  units  assemble  into  a  struc¬ 
ture  resembling  the  crystal  structure  of  NH4CI  promoting 
proton  transfer  and  formation  of  discrete  NH4“^  and  Cl" 
ions.  The  presence  of  water  stabilizes  the  ion-pair  over  the 
hydrogen-bonded  NH3-HCI  complex.  However,  maximum 
stabilization  promoting  the  proton  transfer  is  not  achieved 
until  two  molecules  of  water  are  added  to  the  NH3-HCI 
complex. 

6.2.  Kinetics  of  ammonium  chloride  formation  and 
nucleation 

The  kinetics  of  formation  of  NH4CI  in  the  gas  phase  is 
more  complicated  than  the  simple  homogeneous  reaction  of 
NH3  and  HCl  vapors.  The  particles  formed  grow  by  either 
of  two  mechanisms: 

(1)  By  deposition  of  additional  NH4CI  formed  by  homoge¬ 
neous  reaction  in  the  gas  phase. 

(2)  By  adsorption  of  either  of  the  two  gas-phase  reagents  on 
the  solid  surface  followed  by  gas-solid  reaction  of  the 
adsorbed  species  with  the  other  one  diffusing  in  from 
the  gas  phase  to  the  surface  of  the  particle. 

The  partial  pressure  of  NH3  and  HCl  must  exceed  a 
minimal  critical  value  for  particles  of  ammonium  chloride 
to  be  formed.  Twomey  [165]  found  the  partial  pressure 
of  both  reagents  must  exceed  0.005  mm  of  Hg  for  nu¬ 
cleation  to  start.  Once  the  threshold  partial  pressure  is 
exceeded  the  nucleation  rate  increases  steeply,  increasing 
the  partial  pressure  by  a  factor  of  6  increases  the  nucle¬ 
ation  rate  by  7  orders  of  magnitude.  In  later  experiments, 
Luria  and  Cohen  [166]  found  higher  values  for  the  onset 
of  nucleation,  3.5  x  10^"^  molecules  cm“^  (0.010 mmHg) 
in  a  stationary  system  and  four  to  five  times  higher  values 
still  in  a  flow  system.  It  seems  likely  the  critical  mini¬ 
mal  concentration  for  the  onset  of  nucleation  depends  on 
experimental  conditions  (gas  flow  rates,  cell  geometry, 
temperature). 

Countess  and  Heicklen  [167]  found  the  homogeneous  bi- 
molecular  reaction  rate  for  the  gas-phase  reaction  was  1.9  x 
10"^^  cm^  s"^  The  number  of  particles  produced  increases 
first  linearly  with  time  then  levels  off  as  the  reactants,  NH3 
and  HCl,  are  consumed.  Luria  and  Cohen  [166]  found  the 
same  behavior  proved  growth  of  particles  does  not  occur 
by  particle  coagulation  but  rather  takes  place  by  adsorption 
of  one  of  the  reagents  onto  the  particle  surface  followed  by 
gas-solid  reaction  with  the  other  reagent.  The  processes  of 
nucleation  and  subsequent  particle  growth  are  described  by 


the  following  kinetics  equations: 

Nucleation  :  /?{N} 

=  ^n.2  g-i  [nH3][HC1]2’  (48) 

where  /?{N}  represents  the  particle  production  rate. 

Particle  growth  :  /^{NILtCl} 

=  7  X  10-26  ^.J„6  5-1  [nh3][HC1][N]  (49) 

where  /?{NH4C1}  is  the  rate  of  NH4CI  formation  in  the 
solid  phase  and  [N]  is  the  concentration  of  particles  already 
formed. 

There  is  still  considerable  uncertainty  associated  with  the 
determination  of  the  nucleation  rate  (Eq.  (48)),  experimental 
errors  being  magnified  by  the  derivation  of  the  rate  constant 
from  the  slope  of  a  log-log  plot  of  the  data. 

Under  normal  atmospheric  conditions,  the  partial  pres¬ 
sure  of  both  ammonia  and  HCl  are  below  the  critical  min¬ 
imal  value  for  the  onset  of  nucleation.  It  is  therefore  likely 
that  ammonium  chloride  will  form  on  the  surface  of  exist¬ 
ing  particles  where  ammonia  or  HCl  would  have  adsorbed. 
This  will  certainly  occur  in  the  case  of  a  massive  release  of 
ammonia  during  which  the  high  partial  pressure  of  ammonia 
will  drive  the  reaction. 

63.  Thermodynamics  of  ammonium  chloride  dissociation 

Pio  and  Harrison  [168,169]  conducted  a  critical  review 
and  assessment  of  the  thermodynamics  of  ammonium  chlo¬ 
ride  decomposition.  As  discussed  in  the  preceding  section 
for  ammonium  nitrate,  relative  humidity  and  temperature 
are  the  two  main  determinants  for  the  equilibrium  con¬ 
stant  for  ammonium  chloride  dissociation  according  to 
Eq.  (50). 

NH4C1(s)  NH3  (g)  +  HCl(g)  (50) 

The  point  of  deliquescence  decreases  linearly  from  85% 
RH  at  -10  °C  to  68%  RH  at  70  °C.  Therefore,  under  nor¬ 
mal  atmospheric  conditions,  the  deliquescence  point  will 
range  from  75  to  80%  RH.  Below  this  range,  NH4CI  will 
be  present  as  solid  particles  and  above  as  liquid  droplets 
of  concentrated  aqueous  solution.  Using  the  best  available 
thermodynamic  values  in  Van’s  Hoff  equation,  Pio  and 
Harrison  [168]  derived  Eq.  (51)  for  the  temperature  de¬ 
pendence  of  the  equilibrium  constant  for  solid  ammonium 
chloride: 

ln(Kp)  =  2.2358  ln(7)  -  2.13204  x  10“^  7"* 

+  65.437516  -  8.167  x  lO'^r 
+4.64383  X  lO-’r^  -  1.10475  x  10-‘°^^ 
Kp298K  =  71.56ppb2  (51) 

Values  calculated  using  Eq.  (51)  are  reasonably  close  to 
those  of  the  empirical  Eq.  (52)  derived  by  Wagner  and  New¬ 
man  [170]  from  measurements  of  vapor  pressure  between 
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35  and  85  X: 

ln(Kp)  =  1.812 ln(r)  -  2.138641  x  lO'^r"*  +  68.258553 
-7.4604  X  lO"^?;  Kp298K= 102.3  ppb^  (52) 

As  observed  for  ammonium  nitrate,  above  the  deliques¬ 
cence  point,  the  product  of  partial  pressure  of  ammonia  and 
hydrochloric  acid  above  the  aqueous  droplet  is  much  lower 
than  the  corresponding  product  for  the  solid  particle  below 
the  deliquescence  point.  At  a  given  temperature,  the  partial 
pressure  product  at  95%  RH  is  10-20  times  lower  than  the 
corresponding  product  of  partial  pressure  above  the  solid. 
Increasing  RH  to  98%  further  decreases  the  product  60-100 
times  from  the  value  expected  over  the  solid. 

Measurements  of  atmospheric  concentrations  of  NH3, 
HCl,  and  HNO3  in  England  and  the  Netherlands  [145] 
revealed  that  in  large,  the  measured  concentrations  are  in 
agreement  with  those  predicted  for  equilibrium  chemical 
thermodynamics  of  a  non-interactive  external  mixture  of 
ammonium  chloride  and  nitrate.  Deviations  are  seen  at  low 
temperatures  (<5°C)  and  high  relative  humidity  (>80% 
RH).  Although  no  definitive  explanation  can  be  found,  it 
is  believed  these  deviations  result  from  kinetic  constraints 
preventing  the  system  from  reaching  equilibrium  after  a 
rather  rapid  change  of  atmospheric  conditions. 

6,4,  Kinetics  of  ammonium  chloride  dissociation 

Harrison  et  al.  [146]  measured  the  evaporation  rate  of  am¬ 
monium  nitrate  and  chloride  aerosols  in  an  annular  denuder 
where  vapors  of  NH3,  HCl,  and  HNO3  are  rapidly  and  con¬ 
tinuously  removed.  Rates  were  expressed  in  terms  of  the  rate 
of  decrease  of  particle  diameter.  The  results  are  summarized 
in  Table  14. 

Harrison  et  al.  [146]  reached  the  following  conclusions: 

•  Evaporation  rates  measured  are  independent  of  particle  ra¬ 
dius  and  low  compared  to  those  predicted  by  existing  the¬ 
ories  of  aerosol  evaporation,  implying  that  achievement 
of  equilibrium  is  hindered  by  some  kinetic  constraint.  It 
was  postulated,  the  low  evaporation  rate  of  NH4CI  from 
the  solid  particle  would  be  consistent  with  a  mechanism 
where  the  rate  limiting  step  would  be  the  conversion  of 
NH4CI  into  physically  adsorbed  NH3  and  HCl  at  the  sur¬ 
face  of  the  particle. 

•  Ammonium  chloride  dissociates  faster  than  ammonium 
nitrate,  particularly  from  aqueous  solutions. 


Table  14 

Evaporation  rate  of  ammonium  chloride  and  nitrate  aerosols 


Evaporation  rate  (As  ') 

NH4CI  NH4NO3 

Dry  aerosol 

-1.05  -0.45 

Liquid  aerosol 

-4.52  -0.49 

•  For  liquid  droplets,  the  rate-limiting  step  could  be  the 
diffusion  of  the  solvated  ions  within  the  liquid. 

•  Differences  between  chloride  and  nitrate  evaporation  from 
liquid  droplets  may  be  due  to  differences  in  deliquescence 
points  between  the  two  salts  leading  to  different  solute 
concentrations  at  the  same  RH  and  to  a  difference  in  con¬ 
centration  changes  for  small  RH  changes  in  the  high  rel¬ 
ative  humidity  range. 

7.  Reactions  of  ammonia  with  carbon  dioxide 

Carbon  dioxide  is  one  of  the  main  chemical  components 
of  the  atmosphere  where  concentrations  average  360  ppm. 
Gas-phase  reaction  of  ammonia  and  CO2  is  not  expected  to 
occur  and  has  never  been  reported  under  atmospheric  con¬ 
ditions,  however,  in  aqueous  solutions,  CO2  forms  carbonic 
acid,  a  weak  acid  that  reacts  with  ammonia  to  give  three 
possible  products: 

NH3  +  C02+H20<e>NH4+  +  HC03-  (53) 

NH3  +  NH4+  +  HC03"  2(NH4)+  +  003^"  (54) 

NH3  +  HCO3"  NH2COO-  +  H2O  (55) 

As  indicated  in  Eqs.  (53)-(55),  the  species  formed:  ammo¬ 
nium  bicarbonate,  ammonium  carbonate,  and  ammonium 
carbamate  are  in  equilibrium  with  NH3  and  CO2  vapors. 
Under  normal  atmospheric  conditions,  where  ammonia  con¬ 
centrations  are  in  the  low  ppb  range,  the  contribution  of  these 
reactions  is  essentially  negligible. 

In  the  case  of  a  massive  release  of  ammonia  in  the  atmo¬ 
sphere,  the  high  partial  pressure  of  ammonia  will  lead  to  a 
high  concentration  of  dissolved  ammonia  in  water  droplets  in 
the  atmosphere  as  well  as  in  any  surface  water  in  contact  with 
the  ammonia  cloud.  This  dissolved  ammonia  will  in  turn  re¬ 
act  with  atmospheric  CO2  according  to  Eqs.  (53)-(55). 

The  C02’-NH3-H20  system  is  of  considerable  interest  in 
the  chemical  industry  where  it  is  used,  for  example,  for  the 
scrubbing  of  flue  gases  from  coal  gasification  operations. 
The  most  extensive  set  of  experimental  data  was  developed 
by  van  Krevelen  et  al.  [14]  and  has  been  used  ever  since  to 
test  various  models  describing  the  equilibrium  in  aqueous 
solutions  of  weak  electrolytes  [15,171,172]. 

Van  Krevelen  et  al.  [14]  used  a  simple  Setschenow  cor¬ 
rection  for  non-ideal  solutions  in  which  the  logarithm  of 
the  activity  coefficient  is  assumed  to  be  a  linear  function  of 
the  ionic  strength,  to  calculate  empirical  values  for  the  three 
equilibrium  constants  for  Eqs.  (53)-(55).  Using  these  equi¬ 
librium  constants  and  Henry’s  law  constants,  they  outlined 
an  iterative  method  to  calculate  the  composition  of  the  aque¬ 
ous  solution  and  vapor  pressure  above  it  for  a  wide  range 
of  ammonia  concentrations,  ratio  of  mole-fractions  of  total 
CO2  and  total  ammonia  (R  =  (CO2  +  C03H“  +  CO3”  + 
NH2C00“)/(NH3H-NHJ^ -I-NH2COO”)),  and  temperature. 
The  resulting  graph  for  high  ammonia  concentration  (20  N) 
and  20  °C,  is  shown  in  Fig.  5. 
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Fig.  5.  Equilibrium  of  ammonia  and  carbon  dioxide  over  aqueous  solutions 
[14]. 


The  upper  triangle  in  Fig.  5  shows  the  speciation  for  CO2, 
the  lower  rectangle  shows  the  speciation  for  ammonia.  For 
the  case  of  massive  release  of  ammonia,  R  will  be  close  to 
zero.  The  CO2  in  solution  will  be  present  as  carbonate  and 
bicarbonate  and  very  little  if  any  carbamate  would  be  ex¬ 
pected.  As  long  as  the  partial  pressure  of  ammonia  above 
the  water  surface  is  high,  the  partial  pressure  of  CO2  will 
be  very  small.  The  amount  of  atmospheric  CO2  trapped  in 
solution  as  ammonium  carbonate/bicarbonate  will  be  lim¬ 
ited  by  the  rate  of  gas-phase  diffusion  of  atmospheric  CO2 
toward  the  water  surface.  Given  that  condensed  water  must 
be  present  and  that  the  initial  concentration  of  CO2  is  low, 
360  ppm,  the  overall  impact  of  atmospheric  CO2  on  ammo¬ 
nia  removal  can  be  neglected  for  all  practical  purpose. 


8.  Reactions  of  ammonia  with  atmospheric  oxidants 


2.7  X  10“^^  cm^  molecule”^  s“^  at  room  temperature.  Diau 
et  al.  [173]  critically  reviewed  these  investigations  and  con¬ 
cluded  that  experimental  conditions,  particularly  the  use  of 
Xe  lamps  instead  of  lasers  in  the  flash  photolysis  technique 
to  generate  OH  radicals  and  the  very  low  pressure  used  could 
be  the  sources  of  errors.  In  their  investigations,  they  used  ei¬ 
ther  a  flash  lamp  or  an  UV  laser  to  generate  the  radical  reac¬ 
tant  while  avoiding  the  production  of  NH2  by  photolysis  of 
ammonia  at  200  nm.  A  wide  range  of  pressure  (68-504  Torr) 
and  temperature  (273-433  K)  was  covered  to  ascertain  the 
effect  of  these  variables.  They  demonstrated  that  the  im¬ 
pact  of  photolytically-generated  NH2  was  negligible  under 
their  experimental  conditions  and  that  the  rate  constant 
forEq.  (56)  was  1.47(±0.4)  x  10“^^  cm^molecule“^  s"^ 
in  good  agreement  with  Stuhl  [174]  and  at  least  five 
other  investigators.  The  Arrhenius  equation  for  Eq.  (56)  is 
given  by  Eq.  (57).  The  corresponding  activation  energy  is 
l.SKcalmol"^ 


/:  (cm^  molecule  ^s  ^) 


=  (3.29  ±  1.02)  X  10*“^2exp 


-(922  ±100)1 
T 


J 


(57) 


Based  on  the  careful  selection  of  experimental  conditions 
and  consensus  with  the  largest  number  of  other  independent 
investigations,  the  rate  constant  determined  by  Diau  et  al.  is 
the  most  reliable  value  at  this  point.  This  value  of  the  rate 
constant  would  place  the  average  lifetime  of  ammonia  in  the 
troposphere  at  34  days  if  oxidation  by  OH  radicals  were  the 
only  source  of  removal. 


8.2.  Reaction  with  ozone 


Heicklen  and  co-workers  [175,176]  have  shown  that  am¬ 
monia  reacts  with  ozone  to  give  crystalline  NH4NO3  ac¬ 
cording  to  Eq.  (58). 

2NH3  ±  4O3  ->  NH4NO3  ±  4O2  ±  H2O  (58) 


8.1.  Reaction  with  hydroxyl  radicals 

The  primary  reaction  of  ammonia  with  OH  radicals 
produces  an  amidogen  radical,  NH2,  which  reacts  further 
with  other  oxidative  species  in  the  atmosphere  (O2,  O3, 
OH,  H2O2).  This  primary  reaction  (Eq.  (56))  is  the  rate¬ 
determining  step  in  the  atmospheric  oxidation  of  ammonia 
[173].  The  subsequent  reactions  of  the  amidogen  radical 
have  not  been  investigated  in  detail  but  are  believed  to  result 
in  the  formation  of  NO;c* 

NH3  +  OH  NH2  +  H2O  (56) 

Stuhl  [174]  conducted  the  first  experimental  deter¬ 
mination  of  kinetics  of  Eq.  (56)  and  found  a  rate  con¬ 
stant  of  1.5(±0.4)  X  10“"^^  cm^  molecule”^  s”^  at  298  K. 
Since,  18  separate  determinations  of  this  rate  constant 
have  been  reported  with  values  ranging  from  0.41  to 


The  reaction  proceeds  in  three  major  steps: 

(1)  Induction  period  during  which  no  formation  of  am¬ 
monium  nitrate  is  observed  until  the  products  of  the 
concentrations  of  ammonia  and  nitric  acid  reaches  the 
critical  value  of  5.8  x  10^^  molecules  cm""^.  During 
this  induction  period,  ammonia  is  oxidized  to  nitric 
acid.  The  kinetics  of  this  reaction  is  a  complex  func¬ 
tion  of  the  initial  stoichiometric  ratio  of  ammonia  over 
ozone,  of  the  pressure  of  added  inert  gas,  and  of  the 
surface-to-volume  ratio  of  the  reaction  vessel  used. 

(2)  Particle  nucleation:  Once  the  critical  concentrations  of 
NH3  and  HNO3  are  reached,  particle  formation  occurs 
rapidly.  Analysis  of  the  data  suggests  that  the  initial 
particles  are  crystals  containing  8  monomeric  units 
of  ammonium  nitrate.  The  nucleation  is  completed  in 
less  than  10  min  and  a  particle  density  of  about  10^ 
particles/cm^  is  obtained. 
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(3)  Particle  growth:  The  analysis  of  the  kinetics  of  particle 
growth  indicates  that  growth  by  particle  coagulation 
is  negligible.  The  particles  grow  by  condensation  of 
nitric  acid  vapor  on  the  preformed  particles  followed 
by  reaction  of  the  adsorbed  nitric  acid  with  ammonia 


(Eqs.  (59)  and  (60)). 

HN03+C„  4>C„  HN03 

(59) 

C„HN03+NH3-^C„+i 

(60) 

The  rate  of  particle  growth  is  controlled  by  the  rate  of  con¬ 
densation  of  nitric  acid  on  the  ammonium  nitrate  particles 
(Eq.  (59)).  The  rate  constant  for  reaction  51  was  found  to 
be  1.24  X  10“^cm^min“’^ 

8.3,  Overall  loss  of  ammonia  by  oxidation 

Given  the  very  low  concentrations  of  OH  radicals  (less 
than  1  ppt)  and  ozone  (20-80  ppb)  in  the  troposphere  and 
their  relatively  low  reaction  rates,  the  impact  of  these  oxi¬ 
dants  on  a  massive  local  release  of  ammonia  will  be  negli¬ 
gible.  The  same  holds  true  for  other  oxidative  species  such 
as  H2O2  and  peroxy  radicals  that  are  also  present  at  very 
small  concentrations. 


9.  Interactions  with  airborne  oi^anic  compounds 

Under  atmospheric  conditions,  ammonia  is  not  expected 
to  react  with  most  organic  compounds  found  in  the  atmo¬ 
sphere  either  as  a  result  of  natural  processes  (methane  or  ter- 
penes  from  conifer  trees  for  example)  or  from  anthropogenic 
sources  (hydrocarbons  from  gasoline  for  example).  There  is 
no  indication  in  the  literature  that  atmospheric  reactions  of 
ammonia  with  organic  compounds  occur. 

However,  organic  compounds  may  affect  the  kinetics  of 
reaction  of  ammonia  with  water  and  aqueous  acid  droplets 
in  the  atmosphere.  Rubel  and  Gentry  [177,178]  showed  that 
organic  chemicals  capable  of  forming  a  film  on  the  surface 
of  water  may  have  a  significant  effect  on  the  rate  of  mass 
transfer  and  reaction  of  atmospheric  ammonia  with  water 
and  aqueous  acid  solutions.  They  found  that  covering  an  acid 
droplet  with  a  monolayer  of  solid  hexadecanol  decreases 
the  accommodation  coefficient  of  ammonia  by  a  factor  of 

10.  As  the  monolayer  expands,  i.e.,  the  area  per  monolayer 
molecule  increases,  the  monolayer  experiences  a  phase  tran¬ 
sition  from  solid  to  liquid  state.  The  accommodation  coef¬ 
ficient  decreases  further  until  the  monolayer  collapses  as  it 
contracts  because  of  increased  monolayer  cohesion.  When 
the  monolayer  collapses  a  rapid  increase  in  accommodation 
coefficient  is  observed. 

These  observations  suggest  that  film-forming  organic 
compounds  in  the  atmosphere  are  likely  to  impede  the  dis¬ 
solution  of  ammonia  in  water  and  its  reaction  with  acid 
aerosols.  This  phenomenon  is  not  fully  quantified  and  its 
occurrence  is  difficult  to  predict  as  it  will  depend  on  the  na¬ 


ture  of  the  organic  compounds  present  in  the  atmosphere  at 
a  given  time  and  place.  In  general,  the  prediction  for  ammo¬ 
nia  capture  by  water  and  acid  aerosols  should  be  regarded 
as  upper  limits  only,  the  actual  amount  being  potentially  10 
times  or  more  lower  because  of  the  effect  of  film-forming 
organic  pollutants. 


10«  Modeling  of  the  atmospheric  chemistry  of  ammonia 

Attempts  to  include  the  chemistry  of  ammonia  in  atmo¬ 
spheric  models  fall  into  three  main  categories: 

(1)  Modeling  of  aerosol  formation,  particularly  ammonium 
nitrate  particles. 

(2)  Long-range  modeling  of  ammonia  and  ammonium  salts. 

(3)  Catastrophic  release  of  ammonia  from  industrial  opera¬ 
tions. 

lO.L  Modeling  of  aerosol  formation 

These  models  are  essentially  thermodynamics  analysis 
predicting  the  chemical  composition  and  growth  rate  of  at¬ 
mospheric  aerosols  containing  nitrates,  sulfates,  chlorides, 
ammonium  and  other  species  as  a  function  of  prevailing 
atmospheric  conditions:  temperature,  relative  humidity, 
gas-phase  concentrations  of  the  reactants.  They  are  using 
the  thermodynamic  data  discussed  in  Section  4  for  sulfate,  4 
for  nitrate,  and  5  for  chlorides.  For  example,  the  model  de¬ 
veloped  by  Saxena  and  Seigneur  [179]  accounts  well  for  the 
composition  and  phase  changes  of  aerosol  particles  given 
ambient  relative  humidity  and  gas-phase  concentration  of 
the  reactants. 

Chang  et  al.  [180]  developed  and  used  STEM-II,  a  Eule- 
rian  transport/chemistry  model  that  can  treat  as  many  as  44 
chemical  species,  23  of  which  are  treated  by  the  advection- 
diffiision  equation,  while  the  remaining  species  are  treated 
by  use  of  pseudo-steady  state  thermodynamics.  They  applied 
it  to  predict  the  ammonium  nitrate  aerosol  formation  in  an 
area  of  Japan,  far  and  remote  from  major  sources  of  pol¬ 
lution.  The  model  predicted  accurately  changes  in  aerosol 
composition  as  relative  humidity,  temperature  and  concen¬ 
trations  of  ammonia  and  nitric  acid  changed  during  the  day. 

Campbell  et  al.  [181]  used  EQUISOV  H,  a  model  that 
applied  thermodynamics  to  compute  equilibrium  states  for 
species  distributed  among  a  gas  phase  and  multiple,  inter¬ 
nally  mixed,  size  “bins”  of  an  aerosol  phase.  The  model  was 
generally  successful  in  predicting  the  distribution  of  nitrate, 
chloride,  sulfate,  and  ammonium  among  fine  and  coarse  par¬ 
ticles  collected  over  Tampa  Bay,  FL. 

Kramm  and  Dlugi  [182]  developed  an  elaborate  model 
that  combines  the  chemistry  of  ammonium  nitrate  formation 
and  generally  accepted  description  of  micrometeorological 
phenomena  (transfer  of  momentum,  sensible  heat,  mass)  to 
compute  the  vertical  fluxes  of  nitric  acid,  ammonia,  and  am¬ 
monium  nitrate.  Specifically,  they  use  a  modified  inferen- 
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tial  approach  based  on  constant  flux  approximation  for  con¬ 
served  species  (i.e.,  total  NH3  +  NHj"  and  total  NO3  —  nitric 
acid  +  nitrate)  and  on  solving  the  flux-gradient  equation  by 
using  the  best  available  relationships  for  the  resistance  to 
mass  transfer  of  the  turbulent  region  of  the  atmospheric  sur¬ 
face  layer  and  of  the  underlying  thin  molecular  turbulent  sub 
layer.  The  atmosphere-vegetation  interaction  is  modeled  us¬ 
ing  a  detailed  account  of  the  various  resistances:  stomatal 
resistance  and  cuticular  resistance  of  foliage,  dry  and  wet 
soil  resistances,  etc.  The  model  yields  reasonable  results  for 
the  vertical  distribution  of  chemical  species  when  predict¬ 
ing  expected  changes  caused  by  meteorological  or  chemical 
concentrations.  However,  the  final  test  of  the  model  against 
actual  field  data  is  not  possible  because  of  the  present  limita¬ 
tion  of  available  equipment  to  measure  actual  concentrations 
of  ammonia,  nitric  acid  and  ammonium  nitrate  in  real  time. 

10.2,  Long-range  transport  models 

A  number  of  long-range  transport  models  have  been  de¬ 
veloped  and  tested  in  Europe  to  describe  the  migration  and 
interactions  of  pollutants  across  national  boundaries.  These 
models  generally  include  rate  equations  for  the  reaction  of 
acidic  pollutants  with  ammonia  and  for  the  dry  and  wet 
deposition  of  ammonia  and  anunonium  salts.  Although  the 
general  trends  are  generally  depicted  faithfully,  accuracy  of 
prediction  is  hampered  by  the  lack  of  precise  information  on 
ammonia  emission  and  the  difficulty  of  assessing  whether 
a  particular  area  is  a  net  source  or  sink  for  ammonia  (com¬ 
pensation  point). 

Russell  et  al.  [183]  developed  a  model  describing  the 
formation  and  transport  of  ammonium  nitrate.  Their  La- 
grangian  model  takes  into  account  the  kinetics  of  formation 
of  nitric  acid  from  NO^  and  the  thermodynamics  of  the 
ammonium  nitrate  equilibrium.  The  main  uncertainty  for 
modeling  remains  the  accurate  description  of  ammonia 
emissions  from  natural  and  anthropogenic  sources.  Nev¬ 
ertheless,  the  model’s  predictions  compare  favorably  with 
actual  measurements  made  at  El  Monte,  CA.  It  explains 
satisfactorily  the  peak  of  nitrate  formation  early  in  the 
morning  when  the  HNO3  produced  soon  after  sunrise  reacts 
with  the  ammonia  that  accumulated  during  the  night. 

Long-range  transport  models  developed  to  assess  the  fate 
and  transport  of  ammonia  and  ammonium  ions  in  Europe 
generally  include  selected  average  dry  and  wet  deposition 
velocities  for  ammonia  and  ammonium  ions,  thermodynam¬ 
ics  of  sulfate  and  nitrate  formation  and  approximate  rates 
for  the  oxidation  of  SOx  and  NO^  in  the  atmosphere.  Vari¬ 
ous  approaches  are  used  to  describe  physical  dispersion  of 
the  pollutants:  Lagrangian  [184],  simple  mixing  and  long 
term  averaging  [185],  Gaussian  plume  formulation  for  point 
sources  and  statistical  long  term  transport  [186],  Eulerian 
for  horizontal  transport  and  Gaussian  for  vertical  distribu¬ 
tion  [187,188].  In  general,  long-term  trends  are  described 
satisfactorily  but  the  main  weakness  remains  in  the  lack 
of  good  data  for  local  ammonia  emission.  In  general,  am¬ 


monia  is  emitted  near  the  ground  and  is  transported  only 
short  distances  from  the  point  of  emission.  Ammonium  sul¬ 
fate  and  nitrate  tend  to  be  formed  at  higher  altitude  by  re¬ 
action  of  rising  ammonia  with  acid  pollutants  emitted  at 
higher  elevation  and  transported  often  from  distant  point 
sources.  The  ammonium  salts  are  transported  to  very  distant 
locations. 

Singles  et  al.  [188]  used  a  multi-layer  (33  layers)  model  to 
describe  the  vertical  mixing  of  ammonia  emitted  at  ground 
level  and  selected  values  for  deposition  velocities  to  map 
ammonia/ammonium  concentrations  in  the  British  Isles  at  a 
fine  resolution  level  (5  kmx5  km).  Their  results  were  in  good 
general  agreement  with  measured  values  when  available. 

Apsimon  et  al.  [189]  used  a  Lagrangian  model  simulat¬ 
ing  transport  of  a  column  of  air  divided  in  a  large  num¬ 
ber  of  horizontal  layers  with  finer  resolution  near  ground 
level  to  describe  ammonia  and  ammonium  concentration  in 
one  location  in  England  during  an  anticyclone  episode  in 
June  1983.  The  trajectories  of  air  masses  bringing  pollu¬ 
tants  from  continental  Europe  was  traced  back  and  forma¬ 
tion  of  ammonium  species  from  reactions  of  sulfuric  and 
nitric  acids  produced  by  oxidation  of  SO;c  and  NO;c  with 
ammonia  produced  by  farming  in  the  Netherlands  was  sim¬ 
ulated.  The  results  stress  the  importance  of  limited  mixing 
and  non-homogeneous  concentrations  in  air  masses  and  the 
corresponding  delays  in  chemical  reactions. 

Dentener  and  Crutzen  [190]  built  a  three-dimensional 
global  transport  model  for  ammonia  and  ammonium  salts 
with  a  10°  X  10°  resolution  level.  Transport  is  described  by 
monthly  average  winds  and  eddy  diffusion  parameterization 
based  on  the  standard  deviations  of  wind  from  the  mean.  Re¬ 
moval  of  ammonia  by  reaction  with  acidic  aerosols  (sulfate, 
nitrate),  wet  and  dry  depositions,  and  atmospheric  oxidation 
reactions  are  taken  into  account.  Anthropogenic  emissions 
of  ammonia,  mainly  as  a  result  of  animal  husbandry,  are 
estimated  at  30.4  x  10^  tonnes  N  per  year  exceeding  natu¬ 
ral  sources  by  a  factor  of  two.  The  results  show  reasonable 
agreement  with  measured  values  except  for  Africa  and  Asia 
where  discrepancies  exist. 

10.3,  Catastrophic  releases— modeling  and  simulation 

Release  of  compressed  liquefied  gases  following  the  burst¬ 
ing  of  a  storage  vessel  follows  a  sequence  of  events  that  can 
be  divided  into  a  number  of  distinct  phases: 

•  Initial  flash  expansion,  driven  by  the  stored  internal  en¬ 
ergy,  leading  to  the  formation  of  a  fine  aerosol  and  con¬ 
siderable  turbulence. 

•  Entrainment  of  ambient  air  due  to  the  action  of  this  tur¬ 
bulence,  together  with  evaporation  of  the  droplets,  while 
the  turbulence  decays. 

•  Gravitational  slumping  and  reduction  of  the  rate  of  en¬ 
trainment  of  air  in  the  presence  of  a  stabilizing  density 
gradient. 

•  Turbulent  dispersion  due  to  wind. 
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This  sequence  of  events  is  extraordinarily  complex  and 
has  not  so  far  been  modeled.  In  general,  simulation  and 
modeling  efforts  have  concentrated  on  characterizing  one 
phase  of  the  process  to  eventually  build  a  comprehensive 
model  of  the  event. 

Nielsen  et  al.  [191]  reported  on  a  series  of  controlled 
releases  made  as  part  of  the  European  ‘Tladis  Field 
Experiments”.  The  source  was  a  horizontal  flash  boiling 
jet  with  release  rates  of  0.25-0.5  kg  s"”^.  The  experimen¬ 
tal  site  was  heavily  instrumented  with  five  different  sensor 
types  for  ammonia  concentration  (catalytic  combustion, 
electrochemical,  UV  ionization,  sonic  anemometers  with 
thermocouples,  and  back-scatter  Lidar),  numerous  temper¬ 
ature  measurements,  wind  speed  and  direction,  humidity 
measurements,  and  turbulence  measurements.  These  mea¬ 
surements  gave  a  good  characterization  of  the  release  as  it 
evolves  from  a  two-phase  jet  to  a  plume  with  passive  dis¬ 
persion.  The  plume  was  characterized  by  its  centerline  con¬ 
centration,  height,  and  width  and  the  results  compared  with 
the  prediction  of  several  dispersion  models  and  wind-tunnel 
simulations.  The  concentration  field  was  analyzed  in  a  fixed 
frame  of  reference  as  well  as  in  a  frame  of  reference  mov¬ 
ing  with  the  instantaneous  plume  centerline.  In  general,  the 
moving-frame  statistics  are  in  better  agreement  with  the 
prediction  of  the  dispersion  models  than  the  fixed-frame 
statistics  of  the  field  measurements.  Source  measurements 
of  exit  pressure  and  temperature  allow  estimates  of  release 
enthalpy  and  jet  flow  force.  Evidence  for  heat  input  from  the 
ground  was  found  because  of  the  deviation  of  specific  plume 
enthalpy  near  the  ground  from  expected  adiabatic  mixing  of 
air  and  ammonia.  No  sign  of  plume  lift-off  was  observed. 
Lidar  measurements  demonstrated  the  high  variability  of 
instantaneous  horizontal  plume  profile.  Small-scale  turbu¬ 
lence  within  the  plume  adds  to  plume  meandering  caused 
by  wind  shifts  to  cause  these  concentration  fluctuations.  In 
general,  the  dispersion  models  used  fail  to  describe  plume 
meandering  and  the  effect  of  terrain  topography  (slope  for 
example)  and  nearby  obstacles  (e.g.,  buildings). 

Deaves  et  al.  [192]  developed  a  mathematical  model 
known  as  airborne  concentration  estimate  (ACE)  to  describe 
the  initial  phase  of  the  catastrophic  release  to  complement 
the  standard  dispersion  models.  The  ACE  model  consists 
of  an  explosion  and  of  a  turbulent  growth  sub-model.  The 
explosion  sub-model  takes  as  inputs  the  material  type,  liq¬ 
uid  mass,  storage  and  ambient  temperature  and  pressure 
before  explosion  and  calculates  the  cloud  radius,  tempera¬ 
ture,  masses  of  ammonia  and  air,  and  turbulent  velocity  and 
length  scale  at  the  end  of  the  explosion  phase,  i.e.,  when 
the  main  agent  of  cloud  expansion  changes  from  outward 
momentum  provided  by  the  explosion  to  air  entrainment 
due  to  inherent  violent  turbulence.  The  turbulent  growth 
sub-model  describes  turbulence  decay  from  its  level  at  the 
end  of  the  explosion  phase  to  a  level  comparable  to  the  sum 
of  atmospheric  turbulence  and  turbulence  due  to  gravita¬ 
tional  slumping.  Parameters  in  the  sub-model  are  adjusted 
to  fit  data  obtained  from  small-scale  experiments.  An  energy 


balance  is  established  tracking  the  internal  energy,  kinetic 
energy,  and  work  done  between  the  initial  stage  and  the  end 
of  each  phase.  Liquid  and  gas  mixture  are  assumed  to  be  in 
equilibrium  throughout  the  evaporation  process.  The  model 
was  applied  to  a  range  of  catastrophic  releases,  illustrating 
the  impact  of  the  nature  of  the  compressed  liquid  (chlorine, 
propane,  or  ammonia),  of  the  release  size  (10-100  tonnes), 
and  of  the  release  direction  (upward  or  downward)  on  the 
size  and  composition  of  the  cloud,  aerosol,  and  liquid  pool 
as  well  as  distance  traveled  and  time  before  slumping. 


!!•  Conclusion — recommendations 

The  atmospheric  chemistry  of  ammonia  has  been  exten¬ 
sively  studied  during  the  past  25  years.  Being  the  only  basic 
chemical  naturally  present  in  the  troposphere,  ammonia 
plays  an  important  role  in  the  neutralization  of  acidic  pol¬ 
lutants  such  as  nitrogen  and  sulfur  oxides  and  thus  impacts 
acid  rain  precipitation.  On  the  negative  side,  the  products 
of  these  neutralization  reactions,  ammonium  sulfates  and 
nitrate  principally,  are  the  main  particulate  components 
of  smog  and  thus  significantly  affect  the  opacity  of  the 
atmosphere  and  the  earth  radiation  budget  while  infiuenc- 
ing  precipitation  events  by  serving  as  cloud  condensation 
nuclei. 

Although  the  results  reported  in  the  literature  for  the  reac¬ 
tion  of  ammonia  with  atmospheric  components  are  directly 
applicable  to  the  case  of  a  massive  release  of  ammonia  in 
the  atmosphere,  an  important  qualifier  should  be  kept  in 
mind.  In  the  troposphere,  ammonia  is  present  at  very  low 
concentrations,  ranging  from  1  ppb  in  clean  environments  to 
10-25  ppb  in  polluted  atmosphere.  In  the  event  of  a  massive 
release  of  ammonia,  the  concentration  of  ammonia  will  be 
several  orders  of  magnitude  greater  than  considered  in  the 
reported  studies.  Whereas  the  impact  of  this  increased  con¬ 
centration  can  be  easily  accounted  for  by  using  the  reaction 
rates  and  equilibrium  constants  determined  in  the  literature, 
phenomena  controlled  by  mass  transfer  steps  will  require 
careful  reconsideration.  For  example,  the  mass  transfer  of 
ammonia  from  the  gas  phase  to  the  surface  of  water  con¬ 
trols  the  rate  of  capture  of  ammonia  by  water  under  normal 
atmospheric  conditions,  because  of  the  very  low  ammonia 
concentration  in  the  gas  phase.  Under  conditions  prevailing 
during  a  massive  release  of  ammonia  in  the  atmosphere, 
mass  transfer  from  gas  to  liquid  surface  will  not  be  the 
rate-limiting  step  anymore  because  of  the  high  ammonia 
concentration  in  the  gas  phase.  Under  process  steps  such 
as  the  diffusion  of  dissolved  ammonia  from  the  surface  to 
the  bulk  of  the  liquid  may  become  the  rate-determining 
step. 

Ammonia  readily  dissolves  in  water  where  it  is  partially 
hydrolyzed  to  ammonium  hydroxide  with  an  equilibrium 
constant  =  L774  x  10“^  at  25  °C.  Henry’s  law  governs 
the  vapor  pressure  of  ammonia  over  its  aqueous  solutions 
with  a  constant  that  is  a  function  of  absolute  temperature 
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given  by  the  following  equation: 

/  39 17.507  \ 

IniKn)  =  -8.09694  +  f - - - j  -  0.003147 

The  application  of  Henry’s  law  requires  the  calculation 
of  the  activity  of  dissolved  ammonia.  In  dilute  aqueous 
solutions,  the  activity  of  ammonia  can  be  set  equal  to  its 
molality  without  significant  loss  in  accuracy.  However,  in 
concentrated  solutions  as  may  be  found  in  cloud  droplets  or 
for  aqueous  solutions  containing  significant  concentrations 
of  salts,  such  as  seawater,  the  activity  coefficient  is  best 
computed  by  using  the  Pitzer  model  that  takes  into  account 
the  interaction  of  ammonia  with  dissolved  ions  and  neu¬ 
tral  molecules.  Recommended  values  for  the  corresponding 
interaction  parameters  are  listed  in  Table  5  of  the  paper. 

As  mentioned  earlier,  mass  transfer  phenomena  generally 
control  the  interaction  of  ammonia  and  water  in  the  environ¬ 
ment  and  a  general  quantitative  description  of  the  process  is 
still  the  object  of  some  controversy.  Nevertheless,  to  model 
the  interaction  of  ammonia  with  water  in  the  environment, 
it  is  useful  to  consider  the  following  separate  cases: 

•  In-cloud  scavenging  of  ammonia:  For  the  very  small  drop 
size  characterizing  droplets  within  clouds,  equilibrium  is 
quickly  established  between  gas  phase  and  liquid  phase 
and  the  concentration  of  ammonia  in-cloud  droplets  can 
be  calculated  from  Henry’s  law. 

•  Below-cloud  scavenging:  Raindrop  sizes  are  significantly 
larger  than  in-cloud  droplets  and  their  fall  is  too  short 
for  equilibrium  to  be  established.  Asman  has  developed  a 
model  (Section  2.4.2)  that  calculates  the  amount  of  am¬ 
monia  scavenged  by  rainfall  given  the  temperature  and 
relative  humidity  at  ground  level,  the  concentration  of  am¬ 
monia  below  the  cloud,  the  volume  of  atmosphere  affected 
by  the  rain,  and  the  duration  of  the  rainstorm. 

•  *'Dry  deposition*'  on  ammonia  on  riven  lake,  and  sea  sur¬ 
faces:  The  dry  deposition  of  ammonia  over  natural  bodies 
of  water  is  likely  to  be  a  significant  mechanism  for  am¬ 
monia  removal  from  the  atmosphere.  Available  data  are 
limited  to  the  very  low  ammonia  concentrations  normally 
present  in  the  atmosphere  and  for  which  gas-phase  dif¬ 
fusion  is  the  controlling  process.  For  the  high  gas-phase 
concentrations  of  ammonia  expected  from  massive  release 
incidents,  the  rate-determining  step  is  likely  to  be  very 
different.  Relevant  data  are  not  available  in  the  literature. 
We  recommend  performing  a  series  of  experiments  to  de¬ 
termine  the  deposition  velocity  of  ammonia  over  surfaces 
of  water  and  seawater  in  the  concentration  range  expected 
in  the  case  of  massive  release  of  the  gas. 

The  dry  deposition  of  ammonia  over  vegetation  is  a  com¬ 
plex  phenomenon  that  is  still  poorly  understood  and  difficult 
to  determine  experimentally.  In  the  case  of  a  massive  release 
of  anunonia  in  the  atmosphere,  the  high  concentration  of 
ammonia  in  air  will  become  the  main  driving  force  for  the 
deposition  process.  Aerodynamic  and  gas-phase  diffusion 
resistance  will  become  negligible  and  the  physico-chemical 


phenomena  on  surfaces  will  control  deposition  velocity. 
Given  the  uncertainties  associated  with  the  results  in  the 
literature  and  the  difficulties  of  conducting  meaningful  ad¬ 
ditional  experiments,  we  recommend  the  adoption  of  the 
following  values  for  dry  deposition  velocities: 

•  For  forested  land: 

o  maximum  deposition  velocity  =  60mms”^  under  wet 
conditions  and  at  100%  humidity; 
o  apply  Sutton’s  exponential  decrease  (Eq.  (30)  in  Section 
3.3.4)  to  calculate  the  deposition  velocity  for  lower  rel¬ 
ative  humidity  levels; 

o  for  nighttime  conditions,  the  maximum  deposition  ve¬ 
locity  (100%  humidity)  =  28mms“^ 

•  For  terrains  with  low  aspect  ratio  vegetation:  agricultural 
crops,  pastures,  tundra,  and  wasteland: 

o  maximum  deposition  velocity  =  20mms"^  under  wet 
conditions  and  at  100%  humidity; 
o  influence  of  humidity  level:  use  Sutton’s  exponential 
decrease  for  maximum  deposition  velocity  value; 
o  For  nighttime  conditions  assume  a  maximum  deposi¬ 
tion  velocity  value  of  lOmms^^  at  100%  humidity. 

As  mentioned  earlier,  ammonia  is  the  only  basic  species 
present  in  the  atmosphere  to  neutralize  acidic  pollutants  such 
as  sulfur  and  nitrogen  oxides,  thus  mitigating  the  adverse 
effect  of  acid  rain.  Despite  the  low  concentrations  of  ammo¬ 
nia  (1-25  ppb),  NOx  (0. 1-750  ppb),  SO^c  (1-200  ppb),  HNO3 
(1-25  ppb)  present  in  the  atmosphere,  the  neutralization  re¬ 
actions  are  rapid  and  reach  completion  or  equilibrium  con¬ 
centration  given  sufficient  atmospheric  mixing.  In  the  case 
of  a  massive  release  of  anunonia,  the  very  large  stoichio¬ 
metric  excess  of  ammonia  over  the  acidic  species  will  en¬ 
sure  rapid,  complete  reaction  resulting  in  the  neutralization 
of  all  acidic  pollutants  encountered  as  the  ammonia  cloud 
expands.  The  net  results  will  be: 

•  Rapid  oxidation  of  SO2  to  SO^:  Ammonia  was  found  to 
catalyze  the  oxidation  of  SO2  to  SO3  under  atmospheric 
conditions.  Under  normal  atmospheric  conditions,  the  re¬ 
action  is  complete  in  less  than  5  min. 

•  Neutralization  of  SO3  by  ammonia:  At  very  low  humid¬ 
ity  levels,  anunonia  will  react  rapidly  with  SO3  to  form 
crystalline  particles  of  sulfamic  acid.  At  higher  humidity 
levels,  aerosols  of  ammonium  sulfates  will  be  formed.  For 
the  high  ammonia  concentrations  expected  from  massive 
release  incidents,  neutral  anunonium  sulfate,  (NH4)2S04, 
is  expected  to  be  the  main  product  formed.  At  relative 
humidity  below  80%,  nanoparticles  of  ammonium  sulfate 
will  predominate  with  possibly  some  sulfamic  acid.  Above 
80%  humidity,  the  ammonium  sulfate  particles  will  grow 
by  water  condensation  until  thermodynamic  equilibrium 
is  reached  for  the  system  [(NH4)2S04]  solid-saturated 
aqueous  solution-NH3  gas. 

•  Neutralization  of  nitric  acid  formed  by  atmospheric  ox¬ 
idation  of  NOx:  The  reaction  is  very  rapid  and  leads  to 
the  production  of  ammonium  nitrate  that  tends  to  de- 
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posit  on  pre-existing  ammonium  sulfate  nanoparticles,  re¬ 
sulting  in  the  formation  of  aerosols  of  coarser  particles 
(10-100  |xm).  Under  normal  atmospheric  conditions,  am¬ 
monium  nitrate  is  in  equilibrium  with  the  vapor  pressure 
of  its  constituents,  NH3  and  NO3H.  However,  with  the 
high  gas-phase  concentration  of  ammonia  expected  dur¬ 
ing  massive  release  incidents,  ammonium  nitrate  will  be 
stable.  The  physical  form  and  chemical  composition  of 
the  aerosol  is  a  complex  function  of  the  relative  humidity 
and  of  the  sulfate-to-nitrate  ratio  (Section  5.3). 

•  Neutralization  of  hydrochloric  acid:  This  reaction  is  ex¬ 
pected  to  be  significant  in  marine  environments  where  HCl 
is  formed  by  reaction  of  sea-salt  aerosols  with  strong  acid 
aerosols  such  as  sulfuric  and  nitric  acids.  The  reaction  is 
very  fast  and  results  in  the  formation  of  ammonium  chlo¬ 
ride.  Like  ammonium  nitrate,  ammonium  chloride  con¬ 
tributes  to  the  growth  of  ammonium  sulfate  particles  and 
is  in  equilibrium  with  the  partial  vapor  pressure  of  its  con¬ 
stituents.  While  high  gas-phase  concentrations  of  ammo¬ 
nia  are  prevailing,  ammonium  chloride  is  expected  to  be 
stable. 

•  Neutralization  of  carbon  dioxide:  To  react  with  ammonia, 
CO2  must  be  dissolved  in  water  and  ionized  as  carbon¬ 
ate  or  bicarbonate  ions.  The  amount  of  atmospheric  CO2 
trapped  in  solution,  as  ammonium  carbonate/bicarbonate 
will  be  limited  by  the  gas-phase  diffusion  of  atmospheric 
CO2  toward  the  water  surface.  Given  that  condensed  wa¬ 
ter  must  be  present  and  that  the  initial  concentration  of 
CO2  is  low,  360  ppm,  the  overall  impact  of  atmospheric 
CO2  on  ammonia  removal  can  be  neglected  for  all  prac¬ 
tical  purpose. 

In  summary,  the  acidic  pollutants  in  the  atmosphere  will 
react  very  rapidly  with  ammonia  present  at  high  concentra¬ 
tion  resulting  in  the  formation  of  ammonium  salt  aerosols 
that  will  increase  the  opacity  of  the  atmosphere  locally.  From 
a  modeling  standpoint,  these  acidic  components  of  the  atmo¬ 
sphere  can  be  considered  as  a  sink  for  ammonia,  removing  a 
stoichiometric  amount  of  ammonia  as  the  cloud  of  released 
gas  travels  through  the  atmosphere.  It  should  be  noted  that 
no  data  was  uncovered  for  the  kinetics  of  reaction  of  am¬ 
monia  with  peroxyacetyl  nitrate  a  common  component  of 
polluted  atmospheres.  PAN  is  known  to  hydrolyze  in  alka¬ 
line  solutions,  forming  acetate  and  nitrate  salts.  Given  the 
low  concentrations  of  PAN  in  the  troposphere  (2-70  ppb),  it 
might  be  reasonable  to  assume  that  PAN  constitutes  another 
sink  for  ammonia  through  alkaline  hydrolysis  in  ammonia 
saturated  water  droplets. 

Ammonia  reacts  relatively  slowly  with  OH  radicals  in  the 
atmosphere  where  its  average  lifetime  would  be  34  days 
if  oxidation  by  OH  radicals  were  the  only  source  of  re¬ 
moval.  This  pathway  of  ammonia  removal  can  therefore  be 
neglected.  The  reaction  with  ozone  is  very  rapid  for  high 
concentration  of  ammonia  and  results  in  the  formation  of 
ammonium  nitrate.  Given  the  low  concentration  of  ozone 
in  the  atmosphere  (80  ppb),  ozone  can  be  considered  as  an¬ 


other  small  sink  for  ammonia,  as  suggested  for  atmospheric 
acidic  components. 

Several  models,  incorporating  the  chemistry  of  ammonia 
in  various  degrees,  have  been  developed  to  account  for  the 
long-range  transport  of  ammonium  species  in  the  atmo¬ 
sphere.  In  general,  ammonia  is  only  transported  to  short 
distances  from  the  emission  source  because  of  efficient 
scavenging  by  water  and  vegetation.  Ammonium  salts  on 
the  other  hand  are  transported  over  longer  distances.  Limited 
information  on  the  modeling  and  on  actual  trials  of  catas¬ 
trophic  release  of  ammonia  from  storage  tanks  stress  the  im¬ 
portance  of  physical  and  thermo  dynamical  processes  during 
the  actual  tank  rupture  on  the  volume  and  shape  of  the  dense 
plume  formed.  Additional  research  in  this  area  is  critical  to 
the  development  of  a  realistic  model  of  massive  release  of 
compressed  liquefied  gas  from  ruptured  containment  vessels. 
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